Digital sequence analysis of dna methylation

ABSTRACT

The present invention relates to methods and compositions for determination of and uses of specific methylation patterns indicative of adenoma and carcinoma. In particular, the invention relates to analysis of defined CpG loci that are coordinately methylated in DNAs from cancer and adenoma samples, methods for identifying coordinately methylated loci, and methods of using analysis of coordinately methylated loci in one or more marker regions in the design of assays for adenoma and cancer.

CROSS-REFERENCE TO RELATED APPLICATIONS

The present application is a divisional of U.S. patent application Ser. No. 13/364,978, filed Feb. 2, 2012, which claims the benefit of U.S. Provisional Patent Application Ser. No. 61/438,649, filed Feb. 2, 2011, each of which is incorporated by reference in its entirety.

FIELD OF THE INVENTION

The present invention relates to methods and compositions for determination of and uses of specific methylation patterns indicative of adenoma and carcinoma. In particular, the invention relates to analysis of defined CpG loci that are coordinately methylated in DNAs from cancer and adenoma samples, methods for identifying coordinately methylated loci, and methods of using analysis of coordinately methylated loci in one or more marker regions in the design of assays for adenoma and cancer having improved sensitivity and specificity.

BACKGROUND OF THE INVENTION

In higher order eukaryotes, DNA may be methylated at cytosines located 5′ to guanosine in CpG dinucleotides. This modification has important regulatory effects on gene expression, especially when involving CpG rich areas, known as CpG islands, often found in the promoter regions of genes. While approximately 75% of the CpG sites throughout the human genome are methylated, CpG sites within CpG islands are normally unmethylated, and aberrant methylation of CpG islands has been associated with certain diseases, including cancers. For example, CpG island hypermethylation is associated with transcriptional inactivation of defined tumor suppressor genes in human cancers, e.g., colorectal cancer. Therefore, detection of hypermethylated nucleic acid could indicate susceptibility or onset of various forms of cancers.

Despite indications suggesting a link between CpG island methylator phenotype (CIMP) and cancers (see, e.g., Baylin S B, et al., Adv Cancer Res 1998; 72:141-196 and Jones P A, et al., Nat Rev Genet 2002; 3:415-428), the idea that analysis of methylation status alone could be a useful diagnostic or prognostic tool has been controversial. As discussed by Issa, et al. in an editorial in Gastroenterology 179(3):2005, researchers had mixed results in confirming the link between CI,MP and cancers. Although CIMP was reportedly demonstrated in multiple other malignancies (Shen, I., et al. J Natl Cancer Inst 2002; 94:755-761; Garcia-Manero G, et al., Clin Cancer Res 2002; 8:2217-2224; Toyota M, et al., Blood 2001; 97:2823-2829; Ueki T, et al., Cancer Res 2000; 60:1835-1839; Toyota M, et al., Cancer Res 1999; 59:5438-5442; Strathdee G, et al., Am J Pathol 2001; 158:1121-1127; Abe M, et al., Cancer Res 2005; 65:828-834) and several groups confirmed the original findings using similar markers and technology (Whitehall V L, et al., Cancer Res 2002; 62:6011-6014; van Rijnsoever M, et al., Gut 2002; 51:797-802) other groups were not able to establish such links (Eads C A, et al., Cancer Res 2001; 61:3410-3418; Esteller M, et al., Cancer Res 2000; 60:129-133). As late as 2003, a publication concluded that all methylation events in colorectal cancer were related to aging rather than neoplasia (Yamashita K, et al., Cancer Cell 2003; 4:121-131).

The discrepant results have been attributed in part to the fact that it has been demonstrated that 70% to 80% of aberrant DNA methylation events in colorectal cancer are age-related (Toyota M, et al., Proc Natl Acad Sci USA 1999; 96:8681-8686) and that cancer-linked phenotypes are only clear when these are filtered out. It has also been noted that overly sensitive, non-quantitative methods can overestimate methylation and mask the distinctions between methylation that is associated with cancer and that which is not. Issa states that “methylation events (alone) may not provide the ideal universal cancer marker they were once thought to be because CIMP target genes will not be useful to screen for all colorectal cancers (many false negatives are predicted), and non-CIMP target genes will likely yield a high rate of false-positives because they are also methylated in normal appearing mucosa of older individuals without tumors” (Issa, et al., supra).

One approach to increase the clinical specificity of methylation analyses in cancer detection is to look at multiple marker genes. For example, Zou, et al., examined the methylation status of BMP3, EYA2, ALX4, and vimentin in cancer samples. While methylation levels were significantly higher in both cancer and adenoma than in normal epithelium, for each of the four genes, the sensitivity as determined by receiver operating curves was not significantly improved by combining any or all markers compared with the best single marker. (Zou, et al., Cancer Epidemiol Biomarkers Prev 2007; 16(12):2686).

Zou also looked at neoplasims showing methylation in more than one of the marker genes and found that co-methylation was frequent, with 72% of the cancers and 84% of the adenomas tested showing hypermehtylation in two or more of the genes. Zou reported that methylation of one or more of four (at least one), two or more of four, three or more of four, or four of four of these marker genes was noted in 88%, 72%, 53%, and 41% of 74 cancers and 98%, 84%, 60%, and 39% of 62 adenomas, compared with 24%, 7%, 3%, and 0% of 70 normal epithelia, respectively, demonstrating that although the assay gets progressively more specific as when more genes are included in the comethylation set, the sensitivity declines precipitously.

SUMMARY OF THE INVENTION

The present invention relates to the methods of identifying regions of specific genes and specific regions of genomic nucleic acid useful in the detection of methylation associated with colorectal cancer. Methods comprise, e.g., detecting methylated sequences in, for example, tissue biopsy, stool extract, or other body fluids with improved sensitivity and specificity. In preferred embodiments, the present invention provides methods of methylation analysis comprising identifying methylation loci showing advantageous methylation ratios when methylation in non-normal cells, e.g., cancer or adenoma cells is compared to background methylation in normal cells. In some embodiments, the present invention relates to methods of analyzing methylation at each of several loci in a set of possible methylation sites within a marker sequence, wherein the presence of methylation at all of the loci within the defined set of sites occurs more frequently in cancer and adenoma cells than in normal cells, such that a finding of methylation at all of the loci in the defined subset of loci in a sample is indicative of adenoma or cancer.

In some embodiments, the present invention provides a method of identifying a set of methylated CpG loci in a marker nucleic acid wherein methylation is indicative of adenoma, comprising:

a) determining the methylation status of a defined set of CpG loci in each of a plurality of individual copies of a marker nucleic acid from a plurality of normal samples;

b) determining the methylation status of said defined set of CpG loci in each of a plurality of individual copies of said marker nucleic acid from a plurality of non-normal (e.g., adenoma or cancer) samples to identify a defined subset of CpG loci from within said defined set,

wherein the percentage of individual copies of said marker nucleic acid from said plurality of normal samples that are methylated at all of said CpG loci in said defined subset is less than the percentage of individual copies of said marker nucleic acid from said plurality of non-normal samples that are methylated at all of said CpG loci in said defined subset, and wherein methylation at all of said CpG loci in said defined subset in said marker nucleic acid is indicative of a non-normal state, e.g., adenoma and/or cancer. In certain embodiments, the mean percentage of individual copies of the marker nucleic acid methylated at all loci in said defined set of CpG loci in said plurality of non-normal samples is greater than the mean percentage of individual copies of the marker nucleic acid methylated at all loci in said defined set of CpG loci in the plurality of normal samples. In preferred embodiments, the mean percentage of individual copies of the marker nucleic acid methylated at all loci in said defined set of CpG loci in the plurality of non-normal samples is at least one standard deviation, preferably two standard deviations, more preferably three standard deviations greater than the mean percentage of individual copies of said marker nucleic acid methylated at all loci in said defined set of CpG loci in said plurality of normal samples.

In some embodiments, the defined subset of CpG loci consists of the same loci in the defined set of CpG loci.

Determination of the methylation status of the set of CpG loci may be accomplished by any method known to those of skill in the art. In some embodiments, the method comprises treating DNA from the samples with bisulfite. Bisulfite modification treatment is described, e.g., in U.S. Pat. No. 6,017,704, the entire disclosure of which is incorporated herein by reference. In some embodiments, determining the methylation status of the defined set of CpG loci comprises digital analysis of each of a plurality of CpG loci in a plurality of individual copies of a marker nucleic acid. In some preferred embodiments, digital analysis comprise digital sequencing, and/or digital PCR.

In certain preferred embodiments, non-normal sample comprises an adenoma sample, and in particular preferred embodiments, comprises a colorectal adenoma sample. In some preferred embodiments, a non-normal sampled comprises a cancer sample, and in certain preferred embodiments, comprises a colorectal cancer sample.

The present invention provides methods of detecting cancer or adenoma in a sample, e.g., from a subject. In some embodiments, the present invention provides methods comprising determining the methylation status of each CpG locus in a defined subset of CpG loci in at least one cancer or adenoma marker nucleic acid molecule, wherein methylation at each of the CpG loci in the defined subset of CpG loci in the cancer or adenoma marker nucleic acid molecule is indicative of cancer or adenoma in the sample. In certain preferred embodiments, the defined subset comprises at least three CpG loci, while in some preferred embodiments, the defined subset comprises at least four CpG loci or at least five CpG loci.

In certain embodiments, the determining comprises analysis of the CpG loci in a nucleic acid detection assay configured to determine the methylation status of each of the loci in a single nucleic acid detection assay. In some preferred embodiments, the determining comprises analysis of the CpG loci in a nucleic acid detection assay configured to determine the methylation status of each of said loci in a single reaction mixture. In some embodiments, the nucleic acid detection assay comprises a primer extension assay. In certain preferred embodiments, the nucleic acid detection assay may comprise one or more of a nucleic acid amplification assay, a nucleic acid sequencing assay, a structure-specific cleavage assay, a 5′ nuclease cleavage assay, an invasive cleavage assay and/or a ligation assay.

The methods of the present invention are not limited to the analysis of a single cancer or adenoma marker nucleic acid. For example, in some embodiments, the methylation status of each CpG locus in a defined subset of CpG loci in at least one cancer or adenoma marker nucleic acid molecule comprises analysis of nucleic acid molecules from a plurality of cancer or adenoma markers. In some embodiments, the plurality of cancer or adenoma markers comprises at least three cancer or adenoma markers, while in some embodiments, the plurality comprises at least four cancer or adenoma markers. In some preferred embodiments, the cancer or adenoma markers and nucleic acid molecules are selected from the group comprising Vimentin, BMP3, Septin 9, TFPI2, 2 regions of LRAT, and EYA4 markers and nucleic acid molecules. In some embodiments, the assay methods of the present invention are combined with the analysis of one or more other cancer markers, such as fecal occult blood markers (e.g., hemoglobin, alpha-defensin, calprotectin, al-antitrypsin, albumin, MCM2, transferrin, lactoferrin, and lysozyme).

In certain preferred embodiments of the method described herein, a cancer or adenoma marker nucleic acid molecule comprises a vimentin nucleic acid molecule, and in some particularly preferred embodiments, the defined subset of CpG loci in the vimentin nucleic acid molecule comprises loci 37, 40, and 45.

In certain preferred embodiments of the method described herein, a cancer or adenoma marker nucleic acid molecule comprises a BMP3 nucleic acid molecule, and in some particularly preferred embodiments, the defined subset of CpG loci in the BMP3 nucleic acid molecule comprises loci 34, 53, and 61.

In certain preferred embodiments of the method described herein, a cancer or adenoma marker nucleic acid molecule comprises a Septin9 nucleic acid molecule, and in some particularly preferred embodiments, the defined subset of CpG loci in said Septin9 nucleic acid molecule comprises loci 59, 61, 68, and 70.

In certain preferred embodiments of the method described herein, a cancer or adenoma marker nucleic acid molecule comprises a TFPI2 nucleic acid molecule and in some particularly preferred embodiments, the defined subset of CpG loci in said TFPI2 nucleic acid molecule comprises loci 55, 59, 63, and 67.

In certain preferred embodiments of the method described herein, a cancer or adenoma marker nucleic acid molecule comprises an EYA4 nucleic acid molecule, and in some particularly preferred embodiments, the defined subset of CpG loci in said EYA4 nucleic acid molecule comprises loci 31, 34, 37, and 44.

In certain preferred embodiments of the method described herein, the at least one cancer or adenoma marker or nucleic acid molecule comprises a plurality markers or nucleic acid molecules comprising Vimentin, BMP3, Septin9, and TFPI2 markers or nucleic acid molecules.

The present invention further provides methods of selecting a defined set of CpG loci in a marker nucleic acid wherein methylation is indicative of non-normal status, e.g., adenoma or cancer, comprising a) determining the methylation status of a plurality of CpG loci in each of a plurality of individual copies of a marker nucleic acid from a plurality of normal samples; b) determining the methylation status of the plurality of CpG loci in each of a plurality of individual copies of said marker nucleic acid from a plurality of non-normal (e.g., adenoma or cancer) samples; c) determining methylation ratios for each locus in the plurality of said CpG loci in the marker nucleic acid; and d) selecting a defined set of CpG loci in the marker nucleic acid, wherein the defined set of CpG loci comprises a plurality of CpG loci having advantageous methylation ratios correlating with non-normal status (e.g., adenoma or cancer).

In some embodiments, determining the methylation ratios comprises determining the ratio between the mean methylation at each of the plurality of CpG loci in the normal samples to the mean methylation at each corresponding CpG locus in said plurality of CpG loci in the non-normal samples. In preferred embodiments, the plurality of individual copies of a marker nucleic acid analyzed in the normal and non-normal (e.g., adenoma or cancer) samples comprises at least 10, preferably at least 100, more preferably at least 1000, still more preferably at least 10,000 and still more preferably at least 100,000 copies. The number of copies analyzed is not limited to these whole numbers, but may be any integer above about 10. The number of copies from different sample types, e.g., normal and non-normal need not be equal.

In certain preferred embodiments of the methods of selecting a defined set of CpG loci in a marker nucleic acid described herein, the plurality of normal and non-normal (e.g., adenoma or cancer) samples compared comprises at least 10, preferably at least 25, still more preferably at least 100 samples. The number of samples analyzed in not limited to these whole numbers, but may be any integer above about 10. The number of different samples of the different sample types, e.g., normal and non-normal, need not be equal.

In certain embodiments, the defined set of CpG loci comprises at least three CpG loci, preferably at least four CpG loci, more preferably at least five CpG loci.

Determination of the methylation status of the plurality of CpG loci may be accomplished by any method known to those of skill in the art, including those described in more detail, below. In some embodiments, the method comprises treating DNA from the samples with bisulfite. In some embodiments, determining the methylation status of the defined set of CpG loci comprises digital analysis of each of a plurality of CpG loci in a plurality of individual copies of a marker nucleic acid. In some preferred embodiments, digital analysis comprises digital sequencing, and/or digital PCR. Methods of preparing samples, e.g., stool samples, for analysis are also known in the art. See, e.g., U.S. Pat. Nos. 7,005,266; 6,303,304; 5,741,650; 5,952,178; and 6,268,136, each incorporated herein by reference.

DEFINITIONS

To facilitate an understanding of the present invention, a number of terms and phrases are defined below.

As used herein, the terms “digital sequencing” and “single molecule sequencing” are used interchangeably and refer to determining the nucleotide sequence of individual nucleic acid molecules. Systems for individual molecule sequencing include but are not limited to the 454 FLX™ or 454 TITANIUM™ (Roche), the SOLEXA™/Illumina Genome Analyzer (Illumina), the HELISCOPE™ Single Molecule Sequencer (Helicos Biosciences), and the SOLID™ DNA Sequencer (Life Technologies/Applied Biosystems) instruments), as well as other platforms still under development by companies such as Intelligent Biosystems and Pacific Biosystems.

As used herein, the term “background” as used in reference to methylation of a locus or region refers to methylation observed in a normal cell or sample at a nucleic acid locus or region that is generally unmethylated in normal cells. For example, CpG islands are generally considered unmethylated in normal human cells but methylation is not completely absent in the CpG islands of normal cells.

As used herein, “methylation” or “methylated,” as used in reference to the methylation status of a cytosine, e.g., in a CpG locus, generally refers to the presence or absence of a methyl group at position 5 of the cytosine residue (i.e., whether a particular cytosine is 5-methylcytosine). Methylation may be determined directly, e.g., as evidenced by routine methods for analysis of methylation status of cytosines, e.g., by determining the sensitivity (or lack thereof) of a particular C-residue to conversion to uracil by treatment with bisulfite. For example, a cytosine residue in a sample that is not converted to uracil when the sample is treated with bisulfite in a manner that would be expected to convert that residue if non-methylated (e.g., under conditions in which a majority or all of the non-methylated cytosines in the sample are converted to uracils) may generally be deemed “methylated”.

As used herein, the terms “digital PCR”, “single molecule PCR” and “single molecule amplification” refer to PCR and other nucleic acid amplification methods that are configured to provide amplification product or signal from a single starting molecule. Typically, samples are divided, e.g., by serial dilution or by partition into small enough portions (e.g., in microchambers or in emulsions) such that each portion or dilution has, on average, no more than a single copy of the target nucleic acid. Methods of single molecule PCR are described, e.g., in U.S. Pat. No. 6,143,496, which relates to a method comprising dividing a sample into multiple chambers such that at least one chamber has at least one target, and amplifying the target to determine how many chambers had a target molecule; U.S. Pat. No. 6,391,559; which relates to an assembly for containing and portioning fluid; and U.S. Pat. No. 7,459,315, which relates to a method of dividing a sample into an assembly with sample chambers where the samples are partitioned by surface affinity to the chambers, then sealing the chambers with a curable “displacing fluid.” See also U.S. Pat. No. 6,440,706 and U.S. Pat. No. 6,753,147, and Vogelstein, et al., Proc. Natl. Acad. Sci. USA Vol. 96, pp. 9236-9241, August 1999. See also US 20080254474, describing a combination of digital PCR combined with methylation detection.

As used herein, “sensitivity” as used in reference to a diagnostic assay, e.g., a methylation assay, refers to clinical sensitivity—the proportion of positive samples that give a positive result using a diagnostic assay. Sensitivity is generally calculated as the number of true positives identified by the assay, divided by the sum of the number of true positives and the number of false negatives determined by the assay on known positive samples. Similarly, the term “specificity” refers to the proportion or number of true negatives determined by the assay divided by the sum of the number of true negatives and the number of false positives determined by the assay on known negative sample(s).

As used herein in reference to diagnostic or analysis assays, the term “complementary” refers to different assays that, when used together, provide a more sensitive and/or specific result than can be provided by any one of the different assays used alone.

As used herein, the term “informative” or “informativeness” refers to a quality of a marker or panel of markers, and specifically to the likelihood of finding a marker (or panel of markers) in a positive sample.

The term “sample” as used herein is used in its broadest sense. For example, a sample suspected of containing a human gene or chromosome or sequences associated with a human chromosome may comprise a cell, chromosomes isolated from a cell (e.g., a spread of metaphase chromosomes), genomic DNA (in solution or bound to a solid support such as for Southern blot analysis), RNA (in solution or bound to a solid support such as for Northern blot analysis), cDNA (in solution or bound to a solid support) and the like.

As used herein, the term “CpG island” refers to a genomic DNA region that contains a high percentage of CpG sites relative to the average genomic CpG incidence (per same species, per same individual, or per subpopulation (e.g., strain, ethnic subpopulation, or the like). Various parameters and definitions for CpG islands exist; for example, in some embodiments, CpG islands are defined as having a GC percentage that is greater than 50% and with an observed/expected CpG ratio that is greater than 60% (Gardiner-Garden et al. (1987) J Mol. Biol. 196:261-282; Baylin et al. (2006) Nat. Rev. Cancer 6:107-116; Irizarry et al. (2009) Nat. Genetics 41:178-186; each herein incorporated by reference in its entirety). In some embodiments, CpG islands may have a GC content >55% and observed CpG/expected CpG of 0.65 (Takai et al. (2007) PNAS 99:3740-3745; herein incorporated by reference in its entirety). Various parameters also exist regarding the length of CpG islands. As used herein, CpG islands may be less than 100 bp; 100-200 bp, 200-300 bp, 300-500 bp, 500-750 bp; 750-1000 bp; 1000 or more by in length. In some embodiments, CpG islands show altered methylation patterns relative to controls (e.g., altered methylation in cancer subjects relative to subjects without cancer; tissue-specific altered methylation patterns; altered methylation in stool from subjects with colorectal neoplasia (e.g., colorectal cancer, colorectal adenoma) relative to subjects without colorectal neoplasia). In some embodiments, altered methylation involves hypermethylation. In some embodiments, altered methylation involves hypomethylation.

As used herein, the term “CpG shore” or “CpG island shore” refers to a genomic region external to a CpG island that is or that has potential to have altered methylation patterns (see, e.g., Irizarry et al. (2009) Nat. Genetics 41:178-186; herein incorporated by reference in its entirety). CpG island shores may show altered methylation patterns relative to controls (e.g., altered methylation in cancer subjects relative to subjects without cancer; tissue-specific altered methylation patterns; altered methylation in stool from subjects with colorectal neoplasia (e.g., colorectal cancer, colorectal adenoma) relative to subjects without colorectal neoplasia). In some embodiments, altered methylation involves hypermethylation. In some embodiments, altered methylation involves hypomethylation. CpG island shores may be located in various regions relative to CpG islands (see, e.g., Irizarry et al. (2009) Nat. Genetics 41; 178-186; herein incorporated by reference in its entirety). Accordingly, in some embodiments, CpG island shores are located less than 100 bp; 100-250 bp; 250-500 bp; 500-1000 bp; 1000-1500 bp; 1500-2000 bp; 2000-3000 bp; 3000 bp or more away from a CpG island.

The term “target,” when used in reference to a nucleic acid detection or analysis method, refers to a nucleic acid having a particular sequence of nucleotides to be detected or analyzed, e.g., in a sample suspected of containing the target nucleic acid. In some embodiments, a target is a nucleic acid having a particular sequence for which it is desirable to determine a methylation status. When used in reference to the polymerase chain reaction, “target” generally refers to the region of nucleic acid bounded by the primers used for polymerase chain reaction. Thus, the “target” is sought to be sorted out from other nucleic acid sequences that may be present in a sample. A “segment” is defined as a region of nucleic acid within the target sequence. The term “sample template” refers to nucleic acid originating from a sample that is analyzed for the presence of a target.

As used herein, the term “locus” refers to a particular position, e.g., of a mutation, polymorphism, or a C residue in a CpG dinucleotide, within a defined region or segment of nucleic acid, such as a gene or any other characterized sequence on a chromosome or RNA molecule. A locus is not limited to any particular size or length, and may refer to a portion of a chromosome, a gene, functional genetic element, or a single nucleotide or basepair. As used herein in reference to CpG sites that may be methylated, a locus refers to the C residue in the CpG dinucleotide.

As used herein, the term “methylation ratio” refers to the amount or degree of methylation observed for particular methylation region or locus (e.g., a CpG locus in a marker gene or region) in a plurality of non-normal cells (e.g., cells in a particular disease state, such as cancerous or pre-cancerous cells) compared to the amount or degree of methylation observed for the same region or locus in a plurality of normal cells (e.g., cells that are not in the particular disease state of interest). For example, for a CpG locus showing mean methylation of 8.39889% in a sampling of normal cells and a mean methylation of 74.0771% in a sampling of a plurality of adenoma cells, a methylation ratio may be expressed as the ratio of the means determined for normal cells:adenoma cells, or 0.11348. A methylation ratio need not be expressed in any particular manner or by any particular calculation. By way of example and not limitation, the methylation ratio above may alternatively be expressed, e.g., as 8.39889:74.0771; 8.39889/74.0771; 74.0771:8.39889; as a calculated “fold methylation over background” 8.81987, etc.

As used herein, the term “advantageous methylation ratio” refers to a methylation ratio for a locus at which methylation correlates with a cellular status, e.g., a particular disease state (for example, normal, precancerous, cancerous) that, when compared to other methylation loci correlated with the same disease state, displays a higher percentage methylation in a population of non-normal cells compared to background levels of methylation at the same locus in a population of normal cells. In some instances, certain CpG loci e.g., within a methylation marker sequence, display a much greater signal-to-noise, i.e., degree in methylation compared to background than other loci in the same marker sequence. In other instances, certain disease-associated marker genes or regions display advantageous methylation ratios at some or all loci compared to the methylation ratios observed for some or all loci within another marker sequence.

As used herein, the term “coordinately methylated” is used in reference to methylation loci, e.g., CpG loci in a marker sequence, that exhibit a particular pattern of methylation that correlates with a cellular status, e.g., a particular disease state (for example, normal, precancerous, cancerous). In preferred embodiments, methylation loci that are all methylated in a manner correlated with a disease state may be deemed to be coordinately methylated in cells having that disease state. “Coordinate methylation” is not limited to situations in which all of the coordinated loci are methylated. Any pattern of methylation among a particular set of loci that correlates with a cellular status, including patterns in which all of the coordinate loci are methylated, patterns in which the loci exhibit a reproducible pattern of methylation and non-methylation, and patterns in which none of the loci within the set are methylated are all included within the meaning of “coordinately methylated.”

As used herein, the term “coordinate methylation analysis” is used interchangeably with “multimethylation analysis” and refers to an assay in which the methylation statuses of a plurality of individual methylation loci in a marker sequence, e.g., CpG loci, are determined together. In preferred embodiments, coordinate methylation analysis is performed using a digital/single copy method (e.g., digital sequencing) or an assay method configured to interrogate all of the selected CpG loci on each molecule tested, such that the methylation pattern in each single molecule tested is revealed.

As used herein, the term “defined set” of CpG loci (or other methylation loci) refers to the set of CpG loci in a marker gene or region selected for methylation analysis. A defined set of CpG loci in a marker gene or region may comprise all CpG loci in the gene or region, or it may comprise fewer than all of the loci in that gene or region.

As used herein the term “defined subset” of CpG loci (or other methylation loci) refers to a subset of the defined set of CpG loci in a marker gene or region, the methylation of which has been determined to be indicative of a non-normal state, e.g., adenoma or cancer. For example, in coordinate methylation analysis to determine the presence of colorectal cancer, the methylation status of a defined subset of CpG loci in at least one cancer marker nucleic acid molecule is determined, with simultaneous methylation at all of said CpG loci in the defined subset being indicative of cancer in the sample. A defined subset of CpG loci in a marker gene or region may comprise all CpG loci in the defined set, or it may comprise fewer than all of the loci in the defined set of loci in that gene or region.

As used herein, the term “colorectal cancer” is meant to include the well-accepted medical definition that defines colorectal cancer as a medical condition characterized by cancer of cells of the intestinal tract below the small intestine (e.g., the large intestine (colon), including the cecum, ascending colon, transverse colon, descending colon, and sigmoid colon, and rectum). Additionally, as used herein, the term “colorectal cancer” is meant to further include medical conditions which are characterized by cancer of cells of the duodenum and small intestine (jejunum and ileum).

As used herein, the term “metastasis” is meant to refer to the process in which cancer cells originating in one organ or part of the body relocate to another part of the body and continue to replicate. Metastasized cells subsequently form tumors which may further metastasize. Metastasis thus refers to the spread of cancer from the part of the body where it originally occurs to other parts of the body. As used herein, the term “metastasized colorectal cancer cells” is meant to refer to colorectal cancer cells which have metastasized; colorectal cancer cells localized in a part of the body other than the duodenum, small intestine (jejunum and ileum), large intestine (colon), including the cecum, ascending colon, transverse colon, descending colon, and sigmoid colon, and rectum.

As used herein, “an individual is suspected of being susceptible to metastasized colorectal cancer” is meant to refer to an individual who is at an above-average risk of developing metastasized colorectal cancer. Examples of individuals at a particular risk of developing metastasized colorectal cancer are those whose family medical history indicates above average incidence of colorectal cancer among family members and/or those who have already developed colorectal cancer and have been effectively treated who therefore face a risk of relapse and recurrence. Other factors which may contribute to an above-average risk of developing metastasized colorectal cancer which would thereby lead to the classification of an individual as being suspected of being susceptible to metastasized colorectal cancer may be based upon an individual's specific genetic, medical and/or behavioral background and characteristics.

The term “neoplasm” as used herein refers to any new and abnormal growth of tissue. Thus, a neoplasm can be a premalignant neoplasm or a malignant neoplasm.

The term “neoplasm-specific marker” refers to any biological material that can be used to indicate the presence of a neoplasm. Examples of biological materials include, without limitation, nucleic acids, polypeptides, carbohydrates, fatty acids, cellular components (e.g., cell membranes and mitochondria), and whole cells. In some instances, markers are particular nucleic acid regions, e.g., genes, intragenic regions, etc. Regions of nucleic acid that are markers may be referred to, e.g., as “marker genes,” “marker regions,” “marker sequences,” etc.

The term “colorectal neoplasm-specific marker” refers to any biological material that can be used to indicate the presence of a colorectal neoplasm (e.g., a premalignant colorectal neoplasm; a malignant colorectal neoplasm). Examples of colorectal neoplasm-specific markers include, but are not limited to, exfoliated epithelial markers (e.g., bmp-3, bmp-4, SFRP2, vimentin, septin9, ALX4, EYA4, TFPI2, NDRG4, FOXE1, long DNA, BAT-26, K-ras, APC, melanoma antigen gene, p53, BRAF, and PIK3CA) and fecal occult blood markers (e.g., hemoglobin, alpha-defensin, calprotectin, al-antitrypsin, albumin, MCM2, transferrin, lactoferrin, and lysozyme). See also U.S. Pat. No. 7,485,420; U.S. Pat. No. 7,432,050; U.S. Pat. No. 5,352,775; U.S. Pat. No. 5,648,212; U.S. RE36713; U.S. Pat. No. 5,527,676; U.S. Pat. No. 5,955,263; U.S. Pat. No. 6,090,566; U.S. Pat. No. 6,245,515; U.S. Pat. No. 6,677,312; U.S. Pat. No. 6,800,617; U.S. Pat. No. 7,087,583; and U.S. Pat. No. 7,267,955, each incorporated herein by reference.

Additional markers include but are not limited those in Table 1, below:

TABLE 1 Accession Symbol GeneID Reference NM_000038 APC 324 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_000044 AR 367 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf AB033043 KIAA1217 56243 www.methdb.de/ and/or www.mdanderson.org/ AK055404 KIAA0984 23329 www.methdb.de/ and/or www.mdanderson.org/ AK090480 www.methdb.de/ and/or www.mdanderson.org/ BC041476 www.methdb.de/ and/or www.mdanderson.org/ BX648962 DKFZp686K1684 440034 www.methdb.de/ and/or www.mdanderson.org/ NM_000017 ACADS 35 www.methdb.de/ and/or www.mdanderson.org/ NM_000022 ADA 100; www.methdb.de/ and/or www.mdanderson.org/ 79015 NM_000038 APC 324 www.methdb.de/ and/or www.mdanderson.org/ NM_000038 APC 324 Weber et al. Nature Genetics 37(8), 2005, 853-862 NM_000043 FAS 355; Weber et al. Nature Genetics 37(8), 2005, 853-862 819114 NM_000044 AR 367 www.methdb.de/ and/or www.mdanderson.org/ NM_000044 AR 367 Weber et al. Nature Genetics 37(8), 2005, 853-862 NM_000076 CDKN1C 1028 www.methdb.de/ and/or www.mdanderson.org/ NM_000077 CDKN2A 1029; DNA Methylation And Cancer Therapy, Landes 51198 Bioscience 2005, ed. Moshe Szyf NM_000077 CDKN2A 1029; www.methdb.de/ and/or www.mdanderson.org/ 51198 NM_000077 CDKN2A 1029; Weber et al. Nature Genetics 37(8), 2005, 853-862 51198 NM_000088 COL1A1 1277 www.methdb.de/ and/or www.mdanderson.org/ NM_000095 COMP 1311 www.methdb.de/ and/or www.mdanderson.org/ NM_000104 CYP1B1 1545 www.methdb.de/ and/or www.mdanderson.org/ NM_000115 EDNRB 1910 www.methdb.de/ and/or www.mdanderson.org/ NM_000115 EDNRB 1910 Weber et al. Nature Genetics 37(8), 2005, 853-862 NM_000115 EDNRB 1910 www.methdb.de/ and/or www.mdanderson.org/ NM_000125 ESR1 2099 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_000125 ESR1 2099 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_000125 ESR1 2099 www.methdb.de/ and/or www.mdanderson.org/ NM_000125 ESR1 2099 Weber et al. Nature Genetics 37(8), 2005, 853-862 NM_000182 HADHA 3030 www.methdb.de/ and/or www.mdanderson.org/ NM_000193 SHH 6469 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000249 MLH1 4292 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_000249 MLH1 4292 www.methdb.de/ and/or www.mdanderson.org/ NM_000280 PAX6 5080 www.methdb.de/ and/or www.mdanderson.org/ NM_000280 PAX6 5080 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000280 PAX6 5080 Weber et al. Nature Genetics 37(8), 2005, 853-862 NM_000308 PPGB 5476 www.methdb.de/ and/or www.mdanderson.org/ NM_000314 PTEN 5728 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_000321 RB1 5925 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_000321 RB1 5925 www.methdb.de/ and/or www.mdanderson.org/ NM_000336 SCNN1B 6338 www.methdb.de/ and/or www.mdanderson.org/ NM_000362 TIMP3 7078 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_000362 TIMP3 7078 www.methdb.de/ and/or www.mdanderson.org/ NM_000362 TIMP3 7078 Weber et al. Nature Genetics 37(8), 2005, 853-862 NM_000378 WT1 7490 www.methdb.de/ and/or www.mdanderson.org/ NM_000402 G6PD 2539 www.methdb.de/ and/or www.mdanderson.org/ NM_000438 PAX3 5077 www.methdb.de/ and/or www.mdanderson.org/ NM_000443 ABCB4 5244 www.methdb.de/ and/or www.mdanderson.org/ NM_000453 SLC5A5 6528 www.methdb.de/ and/or www.mdanderson.org/ NM_000453 SLC5A5 6528 www.methdb.de/ and/or www.mdanderson.org/ NM_000475 NR0B1 190 www.methdb.de/ and/or www.mdanderson.org/ NM_000492 CFTR 1080 www.methdb.de/ and/or www.mdanderson.org/ NM_000492 CFTR 1080 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000514 GDNF 2668 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000517 HBA2 3040 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000520 HEXA 3073; Keshet et al. Nature Genetics 38(2), 2006, 149-153 80072 NM_000524 HTR1A 3350 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000551 VHL 7428 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_000551 VHL 7428 www.methdb.de/ and/or www.mdanderson.org/ NM_000551 VHL 7428 Weber et al. Nature Genetics 37(8), 2005, 853-862 NM_000610 CD44 960 Weber et al. Nature Genetics 37(8), 2005, 853-862 NM_000610 CD44 960 www.methdb.de/ and/or www.mdanderson.org/ NM_000612 IGF2 3481; www.methdb.de/ and/or www.mdanderson.org/ 492304 NM_000620 NOS1 4842 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000680 ADRA1A 148 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000717 CA4 762 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000721 CACNA1E 777 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000782 CYP24A1 1591 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000799 EPO 2056 www.methdb.de/ and/or www.mdanderson.org/ NM_000813 GABRB2 2561 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000818 GAD2 2572 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000829 GRIA4 2893 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000830 GRIK1 2897 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000834 GRIN2B 2904 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000843 GRM6 2916 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000852 GSTP1 2950 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_000852 GSTP1 2950 www.methdb.de/ and/or www.mdanderson.org/ NM_000852 GSTP1 2950 Weber et al. Nature Genetics 37(8), 2005, 853-862 NM_000857 GUCY1B3 2983 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000863 HTR1B 3351 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000902 MME 4311 www.methdb.de/ and/or www.mdanderson.org/ NM_000914 OPRM1 4988 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000915 OXT 5020 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000926 PGR 5241 www.methdb.de/ and/or www.mdanderson.org/ NM_000927 ABCB1 5243 www.methdb.de/ and/or www.mdanderson.org/ NM_000959 PTGFR 5737 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_000965 RARB 5915 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_000965 RARB 5915 www.methdb.de/ and/or www.mdanderson.org/ NM_000965 RARB 5915 Weber et al. Nature Genetics 37(8), 2005, 853-862 NM_000997 RPL37 6167 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001001336 CYB5R2 51700 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001001723 TMEM1 7109 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001002295 GATA3 2625 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001003689 L3MBTL2 83746 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001003891 PCQAP 51586 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_001007792 NTRK1 4914 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001008503 OPRM1 4988 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001008504 OPRM1 4988 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001008505 OPRM1 4988 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001009598 RXRG 6258 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001011545 BACH1 571 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001012331 NTRK1 4914 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001013464 LOC401363 401363; www.methdb.de/ and/or www.mdanderson.org/ 441242; 402532 NM_001018084 SLC26A10 65012 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001020658 PUM1 9698; Keshet et al. Nature Genetics 38(2), 2006, 149-153 28997 NM_001024844 CD82 3732 www.methdb.de/ and/or www.mdanderson.org/ NM_001025205 AP2M1 1173 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001025604 ARRDC2 27106 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001033044 GLUL 2752 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001033056 GLUL 2752 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001033518 WIPI2 26100 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001033519 WIPI2 26100 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001033520 WIPI2 26100 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001033952 CALCA 796 www.methdb.de/ and/or www.mdanderson.org/ NM_001036 RYR3 6263 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001043 SLC6A2 6530 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001053 SSTR5 6755 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001059 TACR3 6870 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001063 TF 7018 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001100 ACTA1 58 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001109 ADAM8 101 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001176 ARHGDIG 398 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001186 BACH1 571 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001204 BMPR2 659 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001228 CASP8 841 www.methdb.de/ and/or www.mdanderson.org/ NM_001250 CD40 958 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001257 CDH13 1012 www.methdb.de/ and/or www.mdanderson.org/ NM_001319 CSNK1G2 1455 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001325 CSTF2 1478 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001385 DPYS 1807; Keshet et al. Nature Genetics 38(2), 2006, 149-153 55412 NM_001451 FOXF1 2294 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001454 FOXJ1 2302 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001458 FLNC 2318 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001480 GALR1 2587 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001538 HSF4 3299 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001553 IGFBP7 3490; www.methdb.de/ and/or www.mdanderson.org/ 818325 NM_001572 IRF7 3665 www.methdb.de/ and/or www.mdanderson.org/ NM_001628 AKR1B1 231 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001635 AMPH 273 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001651 AQP5 362 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001718 BMP6 654 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_001753 CAV1 857 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_001753 CAV1 857 www.methdb.de/ and/or www.mdanderson.org/ NM_001768 CD8A 925 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001801 CDO1 1036 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001851 COL9A1 1297 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001883 CRHR2 1395 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001884 HAPLN1 1404 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001927 DES 1674 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001954 DDR1 780 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_001958 EEF1A2 1917 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001972 ELA2 1991 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001975 ENO2 2026 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_001989 EVX1 2128 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_002007 FGF4 2249 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_002012 FHIT 2272; www.methdb.de/ and/or www.mdanderson.org/ 246734 NM_002024 FMR1 2332 www.methdb.de/ and/or www.mdanderson.org/ NM_002065 GLUL 2752 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_002110 HCK 3055 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_002127 HLA-G 3135 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_002148 HOXD10 3236 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_002155 HSPA6 3310 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_002191 INHA 3623 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_002212 ITGB4BP 3692 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_002221 ITPKB 3707 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_002235 KCNA6 3742 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_002253 KDR 3791 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_002344 LTK 4058 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_002412 MGMT 4255 www.methdb.de/ and/or www.mdanderson.org/ NM_002457 MUC2 4583 www.methdb.de/ and/or www.mdanderson.org/ NM_002478 MYOD1 4654 www.methdb.de/ and/or www.mdanderson.org/ NM_002529 NTRK1 4914 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_002588 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_002658 PLAU 5328; www.methdb.de/ and/or www.mdanderson.org/ 414236 NM_002700 POU4F3 5459 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_002807 PSMD1 5707; Keshet et al. Nature Genetics 38(2), 2006, 149-153 7410 NM_002848 PTPRO 5800 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_002873 RAD17 5884 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_002923 RGS2 5997 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003027 SH3GL3 6457 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003088 FSCN1 6624 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003097 SNRPN 6638; www.methdb.de/ and/or www.mdanderson.org/ 8926; 145624; 8123; 63968; 3653 NM_003149 STAC 6769 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003204 NFE2L1 4779 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003219 TERT 7015 www.methdb.de/ and/or www.mdanderson.org/ NM_003238 TGFB2 7042 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003238 TGFB2 7042 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003246 THBS1 7057 www.methdb.de/ and/or www.mdanderson.org/ NM_003274 TMEM1 7109 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003277 CLDN5 7122 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003295 TPT1 7178; Keshet et al. Nature Genetics 38(2), 2006, 149-153 51447; 2982 NM_003300 TRAF3 7187 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003391 WNT2 7472 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003392 WNT5A 7474 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003408 ZFP37 7539; Keshet et al. Nature Genetics 38(2), 2006, 149-153 7551 NM_003417 ZNF264 9422 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003426 ZNF74 7625 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003435 ZNF134 7693 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003451 ZNF177 7730 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003474 ADAM12 8038 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003508 FZD9 8326 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003539 HIST1H4D 8360 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003540 HIST1H4F 8361 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003541 HIST1H4K 8362 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003546 HIST1H4L 8368 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003666 BLZF1 8548 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003735 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_003736 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_003775 EDG6 8698 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003777 DNAH11 8701; Keshet et al. Nature Genetics 38(2), 2006, 149-153 9026 NM_003806 HRK 8739 www.methdb.de/ and/or www.mdanderson.org/ NM_003823 TNFRSF6B 8771; Keshet et al. Nature Genetics 38(2), 2006, 149-153 51750; 10139 NM_003888 ALDH1A2 8854 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003914 CCNA1 8900 www.methdb.de/ and/or www.mdanderson.org/ NM_003923 FOXH1 8928 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003984 SLC13A2 9058 www.methdb.de/ and/or www.mdanderson.org/ NM_003991 EDNRB 1910 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_003999 OSMR 9180 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004004 GJB2 2706 www.methdb.de/ and/or www.mdanderson.org/ NM_004064 CDKN1B 1027 www.methdb.de/ and/or www.mdanderson.org/ NM_004068 AP2M1 1173 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004102 FABP3 2170 www.methdb.de/ and/or www.mdanderson.org/ NM_004113 FGF12 2257 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004122 GHSR 2693 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004135 IDH3G 3421 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004181 UCHL1 7345 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004230 EDG5 9294 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004248 PRLHR 2834 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004267 CHST2 9435 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004291 CART 9607 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004297 GNA14 9630 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004327 BCR 613; www.methdb.de/ and/or www.mdanderson.org/ 26226 NM_004360 CDH1 999 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_004360 CDH1 999 www.methdb.de/ and/or www.mdanderson.org/ NM_004360 CDH1 999 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_004378 CRABP1 1381 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004385 CSPG2 1462 www.methdb.de/ and/or www.mdanderson.org/ NM_004387 NKX2-5 1482 www.methdb.de/ and/or www.mdanderson.org/ NM_004394 DAP 1611 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004411 DYNC1I1 1780 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004441 EPHB1 2047 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004464 FGF5 2250 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004477 FRG1 2483 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004480 FUT8 2530 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004484 GPC3 2719 www.methdb.de/ and/or www.mdanderson.org/ NM_004525 LRP2 4036 www.methdb.de/ and/or www.mdanderson.org/ NM_004530 MMP2 4313 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004612 TGFBR1 7046 www.methdb.de/ and/or www.mdanderson.org/ NM_004621 TRPC6 7225 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004714 DYRK1B 9149 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004737 LARGE 9215 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004787 SLIT2 9353 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004817 TJP2 9414 www.methdb.de/ and/or www.mdanderson.org/ NM_004865 TBPL1 9519 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004887 CXCL14 9547 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004929 CALB1 793 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004936 CDKN2B 1030 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_004936 CDKN2B 1030 www.methdb.de/ and/or www.mdanderson.org/ NM_004938 DAPK1 1612 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_004975 KCNB1 3745 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004976 KCNC1 3746 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_004988 MAGEA1 4100 www.methdb.de/ and/or www.mdanderson.org/ NM_005032 PLS3 5358 www.methdb.de/ and/or www.mdanderson.org/ NM_005048 PTHR2 5746 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005073 SLC15A1 6564 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005100 AKAP12 9590 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005117 FGF19 9965 www.methdb.de/ and/or www.mdanderson.org/ NM_005117 FGF19 9965 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005157 ABL1 25 www.methdb.de/ and/or www.mdanderson.org/ NM_005159 ACTC 70 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005181 CA3 761 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005233 EPHA3 2042 www.methdb.de/ and/or www.mdanderson.org/ NM_005284 GPR6 2830 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005285 NPBWR1 2831 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005286 NPBWR2 2832 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005346 HSPA1B 3304 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005382 NEF3 4741 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005386 NNAT 4826 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005398 PPP1R3C 5507 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005427 TP73 7161 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_005437 NCOA4 8031 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_005523 HOXA11 3207 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005556 KRT7 3855 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005584 MAB21L1 4081 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005638 SYBL1 6845 www.methdb.de/ and/or www.mdanderson.org/ NM_005668 ST8SIA4 7903 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005806 OLIG2 10215 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005825 RASGRP2 10235 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_005946 MT1A 4489 www.methdb.de/ and/or www.mdanderson.org/ NM_005959 MTNR1B 4544 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006000 TUBA1 7277; www.methdb.de/ and/or www.mdanderson.org/ 84854 NM_006019 TCIRG1 10312 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006041 HS3ST3B1 9953; Keshet et al. Nature Genetics 38(2), 2006, 149-153 84815 NM_006043 HS3ST2 9956 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006053 TCIRG1 10312 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006142 SFN 2810 www.methdb.de/ and/or www.mdanderson.org/ NM_006158 NEFL 4747 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006158 NEFL 4747 www.methdb.de/ and/or www.mdanderson.org/ NM_006161 NEUROG1 4762 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006194 PAX9 5083 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006211 PENK 5179 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006306 SMC1L1 8243 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006307 SRPX 8406 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006463 STAMBP 10617 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006483 DYRK1B 9149 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006484 DYRK1B 9149 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006497 HIC1 3090 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_006497 HIC1 3090 www.methdb.de/ and/or www.mdanderson.org/ NM_006539 CACNG3 10368 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006587 CORIN 10699 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006614 CHL1 10752 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006735 HOXA2 3199 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006765 TUSC3 7991 www.methdb.de/ and/or www.mdanderson.org/ NM_006788 RALBP1 10928 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006874 ELF2 1998; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26472 NM_006898 HOXD3 3232 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_006917 RXRG 6258 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_007117 TRH 7200 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_007181 MAP4K1 11184 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_007182 RASSF1 11186 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_007182 RASSF1 11186 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_007182 RASSF1 11186 www.methdb.de/ and/or www.mdanderson.org/ NM_007197 FZD10 11211 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_007294 BRCA1 672 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_007294 BRCA1 672 www.methdb.de/ and/or www.mdanderson.org/ NM_007332 TRPA1 8989 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_007345 ZNF236 7776 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_007361 NID2 22795 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_012200 B3GAT3 26229 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_012202 GNG3 2785 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_012261 C20orf103 24141 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_012295 CABIN1 23523 www.methdb.de/ and/or www.mdanderson.org/ NM_012295 CABIN1 23523 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_012301 MAGI2 9863 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_012309 SHANK2 22941 www.methdb.de/ and/or www.mdanderson.org/ NM_012309 SHANK2 22941 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_012399 PITPNB 23760 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_012444 SPO11 23626 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_012458 TIMM13 26517 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_013250 ZNF215 7762 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_013291 CPSF1 29894 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_013381 TRHDE 29953 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_013435 RAX 30062 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_013942 PAX3 5077 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014020 LR8 28959 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014080 DUOX2 50506 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014155 BTBD15 29068 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014228 SLC6A7 6534 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014234 HSD17B8 7923 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014325 CORO1C 23603 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014379 KCNV1 27012 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014386 PKD2L2 27039 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014459 PCDH17 27253; Keshet et al. Nature Genetics 38(2), 2006, 149-153 144997 NM_014468 VENTX 27287 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014522 PCDH11X 27328; Keshet et al. Nature Genetics 38(2), 2006, 149-153 83259 NM_014574 STRN3 29966 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014587 SOX8 30812 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014588 VSX1 30813 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014618 DBC1 1620 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014618 DBC1 1620 DNA Methylation And Cancer Therapy, Landes Bioscience 2005, ed. Moshe Szyf NM_014631 SH3PXD2A 9644 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014653 KIAA0789 9671 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014676 PUM1 9698; Keshet et al. Nature Genetics 38(2), 2006, 149-153 28997 NM_014710 GPRASP1 9737 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014724 ZNF96 9753 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014786 ARHGEF17 9828 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014817 KIAA0644 9865 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_014979 SV2C 22987 www.methdb.de/ and/or www.mdanderson.org/ NM_015002 FBXO21 23014 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_015094 HIC2 23119 www.methdb.de/ and/or www.mdanderson.org/ NM_015101 GLT25D2 23127 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_015163 TRIM9 114088 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_015472 WWTR1 25937 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_015507 EGFL6 25975 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_015610 WIPI2 26100 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_015641 TES 26136 www.methdb.de/ and/or www.mdanderson.org/ NM_015683 ARRDC2 27106 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_015722 DRD1IP 50632 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_015920 RPS27L 51065 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016003 WIPI2 26100 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016135 ETV7 51513 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016157 TRO 7216 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016162 ING4 51147 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016179 TRPC4 7223 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016192 TMEFF2 23671 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016192 TMEFF2 23671 www.methdb.de/ and/or www.mdanderson.org/ NM_016223 PACSIN3 29763 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016229 CYB5R2 51700 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016301 ATPBD1C 51184 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016442 ARTS-1 51752 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016521 TFDP3 51270 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016535 ZNF581 51545 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016538 SIRT7 51547 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016540 GPR83 10888 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016552 ANKMY1 51281 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016568 RLN3R1 51289 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016605 FAM53C 51307; Keshet et al. Nature Genetics 38(2), 2006, 149-153 995 NM_016931 NOX4 50507 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016950 SPOCK3 50859 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_016954 TBX22 50945 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_017514 PLXNA3 55558 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_017649 CNNM2 54805 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_017729 EPS8L1 54869 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_017798 YTHDF1 54915 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_017844 ANKMY1 51281 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_017847 C1orf27 54953 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_018061 PRPF38B 55119 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_018074 FLJ10374 55702 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_018129 PNPO 55163 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_018135 MRPS18A 55168 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_018197 ZFP64 55734 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_018310 BRF2 55290 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_018354 C20orf46 55321 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_018401 STK32B 55351 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_018431 DOK5 55816 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_018898 PCDHA6 56142; Keshet et al. Nature Genetics 38(2), 2006, 149-153 56145; 56134; 56147; 56146; 56139 NM_018899 PCDHA6 56142; Keshet et al. Nature Genetics 38(2), 2006, 149-153 56145; 56134; 56147; 56146; 56139 NM_018901 PCDHA6 56142; Keshet et al. Nature Genetics 38(2), 2006, 149-153 56145; 56134; 56147; 56146; 56139 NM_018906 PCDHA6 56142; Keshet et al. Nature Genetics 38(2), 2006, 149-153 56145; 56134; 56147; 56146; 56139 NM_018911 PCDHA6 56142; Keshet et al. Nature Genetics 38(2), 2006, 149-153 56145; 56134; 56147; 56146; 56139 NM_018920 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_018925 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_018926 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_018927 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_018928 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_018950 HLA-F 3134 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_018976 SLC38A2 54407 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_018997 MRPS21 54460 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_019043 APBB1IP 54518 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_019102 HOXA5 3202 www.methdb.de/ and/or www.mdanderson.org/ NM_020166 MCCC1 56922 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_020201 NT5M 56953 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_020208 SLC6A20 54716 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_020226 PRDM8 56978 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_020230 PPAN 56342 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_020348 CNNM1 26507 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_020469 ABO 28 www.methdb.de/ and/or www.mdanderson.org/ NM_020549 CHAT 1103 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_020630 RET 5979 www.methdb.de/ and/or www.mdanderson.org/ NM_020650 RCN3 57333 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_020657 ZNF304 57343 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_020660 CX36 57369 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_020685 C3orf14 57415 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_020815 PCDH10 57575 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_020873 LRRN1 57633 www.methdb.de/ and/or www.mdanderson.org/ NM_020984 CHAT 1103 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_020985 CHAT 1103 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_020986 CHAT 1103 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_020999 NEUROG3 50674 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_021032 FGF12 2257 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_021101 CLDN1 9076 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_021179 C1orf114 57821 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_021193 HOXD12 3238 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_021216 ZNF71 58491 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_021257 NGB 58157 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_021614 KCNN2 3781 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_021911 GABRB2 2561 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_021926 ALX4 60529 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_021956 GRIK2 2898 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_022076 DUSP21 63904 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_022088 ZFP64 55734 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_022169 ABCG4 64137 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_022405 SLC6A20 54716 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_022443 MLF1 4291 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_022468 MMP25 64386; Keshet et al. Nature Genetics 38(2), 2006, 149-153 4328 NM_022469 GREM2 64388 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_022718 MMP25 64386; Keshet et al. Nature Genetics 38(2), 2006, 149-153 4328 NM_022750 PARP12 64761 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_023926 ZNF447 65982 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_024012 HTR5A 3361 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_024046 CAMKV 79012 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_024101 MLPH 79083 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_024306 FA2H 79152 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_024409 NPPC 4880 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_024504 PRDM14 63978 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_024593 EFCAB1 79645 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_024600 C16orf30 79652 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_024826 ASAP 79884 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_024882 C6orf155 79940 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_024893 C20orf39 79953 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_024944 CHODL 140578 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_025019 TUBA4 80086 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_025058 TRIM46 80128 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_025061 LRRC8E 80131 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_025087 FLJ21511 80157 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_025197 CDK5RAP3 80279 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_025204 RP3-402G11.12 80305 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_025208 PDGFD 80310; Keshet et al. Nature Genetics 38(2), 2006, 149-153 414301 NM_025218 ULBP1 80329 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_025263 PRR3 80742 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_030577 MGC10993 80775 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_030667 PTPRO 5800 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_030668 PTPRO 5800 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_030669 PTPRO 5800 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_030670 PTPRO 5800 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_030671 PTPRO 5800 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_030760 EDG8 53637 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_030806 C1orf21 81563; Keshet et al. Nature Genetics 38(2), 2006, 149-153 116492 NM_030920 ANP32E 81611 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_031277 RNF17 56163 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_031283 TCF7L1 83439 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_031424 C20orf55 83541 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_031466 NIBP 83696 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_031488 L3MBTL2 83746 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_031497 PCDHA6 56142; Keshet et al. Nature Genetics 38(2), 2006, 149-153 56145; 56134; 56147; 56146; 56139 NM_031856 PCDHA6 56142; Keshet et al. Nature Genetics 38(2), 2006, 149-153 56145; 56134; 56147; 56146; 56139 NM_031859 PCDHA6 56142; Keshet et al. Nature Genetics 38(2), 2006, 149-153 56145; 56134; 56147; 56146; 56139 NM_031860 PCDHA6 56142; Keshet et al. Nature Genetics 38(2), 2006, 149-153 56145; 56134; 56147; 56146; 56139 NM_031882 PCDHA6 56142; Keshet et al. Nature Genetics 38(2), 2006, 149-153 56145; 56134; 56147; 56146; 56139 NM_031883 PCDHA6 56142; Keshet et al. Nature Genetics 38(2), 2006, 149-153 56145; 56134; 56147; 56146; 56139 NM_031901 MRPS21 54460 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_031912 SYT15 83849 www.methdb.de/ and/or www.mdanderson.org/ NM_031922 REPS1 85021 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_031934 RAB34 83871 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_031994 RNF17 56163 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032034 SLC4A11 83959 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032087 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_032094 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_032098 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_032099 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_032100 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_032101 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_032109 OTP 23440 www.methdb.de/ and/or www.mdanderson.org/ NM_032134 QRICH2 84074 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032140 C16orf48 84080 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032192 PPP1R1B 84152 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032256 TMEM117 84216 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032303 HSDL2 84263 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032391 PRAC 84366 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032402 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_032403 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_032406 PCDHGC3 5098; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26025; 56108; 56112; 9708; 56109 NM_032411 ECRG4 84417 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032412 ORF1-FL49 84418 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032603 LOXL3 84695 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032625 C7orf13 129790 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032803 SLC7A3 84889 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032825 ZNF382 84911 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032838 ZNF566 84924 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032883 C20orf100 84969 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032918 RERG 85004 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032945 TNFRSF6B 8771; Keshet et al. Nature Genetics 38(2), 2006, 149-153 51750; 10139 NM_032961 PCDH10 57575 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_032967 PCDH11X 27328; Keshet et al. Nature Genetics 38(2), 2006, 149-153 83259 NM_032968 PCDH11X 27328; Keshet et al. Nature Genetics 38(2), 2006, 149-153 83259 NM_032969 PCDH11X 27328; Keshet et al. Nature Genetics 38(2), 2006, 149-153 83259 NM_033126 PSKH2 85481 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_033135 PDGFD 80310; Keshet et al. Nature Genetics 38(2), 2006, 149-153 414301 NM_033143 FGF5 2250 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_033224 PURB 5814 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_033302 ADRA1A 148 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_033303 ADRA1A 148 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_033304 ADRA1A 148 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_033445 HIST3H2A 92815 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_033624 FBXO21 23014 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_052902 STK11IP 114790 www.methdb.de/ and/or www.mdanderson.org/ NM_052954 CYYR1 116159 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_052961 SLC26A8 116369 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_052978 TRIM9 114088 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_054021 GPR101 83550 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_054108 HRASLS5 117245 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_058165 MOGAT1 116255 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_078485 COL9A1 1297 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_080552 SLC32A1 140679 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_080617 CBLN4 140689 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_080671 KCNE4 23704 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_080742 B3GAT2 135152 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_130773 CNTNAP5 129684 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_130900 RAET1L 154064 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_133180 EPS8L1 54869 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_133266 SHANK2 22941 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_133338 RAD17 5884 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_133339 RAD17 5884 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_133340 RAD17 5884 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_133341 RAD17 5884 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_133342 RAD17 5884 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_133343 RAD17 5884 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_133344 RAD17 5884 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_133489 SLC26A10 65012 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_133493 CD109 135228 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_133642 LARGE 9215 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_138290 RPIB9 154661 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_138718 SLC26A8 116369 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_138996 CNTNAP5 129684 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_139204 EPS8L1 54869 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_139316 AMPH 273 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_144497 AKAP12 9590 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_144725 FLJ25439 153657 Weber et al. Nature Genetics 37(8), 2005, 853-862 NM_145725 TRAF3 7187 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_145726 TRAF3 7187 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_152562 CDCA2 157313 Weber et al. Nature Genetics 37(8), 2005, 853-862 NM_152854 CD40 958 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_153819 RASGRP2 10235 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_170696 ALDH1A2 8854 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_170697 ALDH1A2 8854 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_170775 KCNN2 3781 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_171827 CD8A 925 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_172337 OTX2 5015 www.methdb.de/ and/or www.mdanderson.org/ NM_173479 LOC126248 126248 www.methdb.de/ and/or www.mdanderson.org/ NM_174869 IDH3G 3421 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_175052 ST8SIA4 7903 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_175611 GRIK1 2897 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_175709 CBX7 23492 www.methdb.de/ and/or www.mdanderson.org/ NM_175768 GRIK2 2898 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_176095 CDK5RAP3 80279 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_176096 CDK5RAP3 80279 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_177555 TRO 7216 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_177556 TRO 7216 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_177557 TRO 7216 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_177558 TRO 7216 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_177959 DOK5 55816 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_178154 FUT8 2530 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_178155 FUT8 2530 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_178156 FUT8 2530 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_178157 FUT8 2530 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_181457 PAX3 5077 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_181458 PAX3 5077 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_181459 PAX3 5077 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_181460 PAX3 5077 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_181461 PAX3 5077 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_181466 ITGB4BP 3692 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_181467 ITGB4BP 3692 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_181468 ITGB4BP 3692 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_181469 ITGB4BP 3692 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_181505 PPP1R1B 84152 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_181657 LTB4R 1241 www.methdb.de/ and/or www.mdanderson.org/ NM_181689 NNAT 4826 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_182609 ZNF677 342926 Weber et al. Nature Genetics 37(8), 2005, 853-862 NM_198265 SPO11 23626 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_198287 ING4 51147 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_198407 GHSR 2693 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_198570 UNQ739 375567 www.methdb.de/ and/or www.mdanderson.org/ NM_198849 LOC283514 283514 Weber et al. Nature Genetics 37(8), 2005, 853-862 NM_199051 FAM5C 339479 www.methdb.de/ and/or www.mdanderson.org/ NM_199076 CNNM2 54805 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_199077 CNNM2 54805 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_199231 GDNF 2668 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_199234 GDNF 2668 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_199425 VSX1 30813 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_199426 ZFP64 55734 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_199427 ZFP64 55734 www.methdb.de/ and/or www.mdanderson.org/ NM_199427 ZFP64 55734 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_201647 STAMBP 10617 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_201999 ELF2 1998; Keshet et al. Nature Genetics 38(2), 2006, 149-153 26472 NM_206827 RASL11A 387496 Weber et al. Nature Genetics 37(8), 2005, 853-862 NM_206866 BACH1 571 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_206961 LTK 4058 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_207121 C20orf55 83541 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NM_213622 STAMBP 10617 Keshet et al. Nature Genetics 38(2), 2006, 149-153 NP_536846 www.methdb.de/ and/or www.mdanderson.org/ NR_002196 www.methdb.de/ and/or www.mdanderson.org/

See also Ilana Keshet, et al., Nature Genetics 38, 149-153 (1 Feb. 2006) and Gerd P Pfeifer, et al., Expert Opinion on Medical Diagnostics, September 2007, Vol. 1, No. 1, Pages 99-108, each of which is incorporated herein by reference.

As used herein, the term “adenoma” refers to a benign tumor of glandular origin. Although these growths are benign, over time they may progress to become malignant. As used herein the term “colorectal adenoma” refers to a benign colorectal tumor in which the cells form recognizable glandular structures or in which the cells are clearly derived from glandular epithelium.

The term “amplifying” or “amplification” in the context of nucleic acids refers to the production of multiple copies of a polynucleotide, or a portion of the polynucleotide, typically starting from a small amount of the polynucleotide (e.g., a single polynucleotide molecule), where the amplification products or amplicons are generally detectable. Amplification of polynucleotides encompasses a variety of chemical and enzymatic processes. The generation of multiple DNA copies from one or a few copies of a target or template DNA molecule during a polymerase chain reaction (PCR) or a ligase chain reaction (LCR; see, e.g., U.S. Pat. No. 5,494,810; herein incorporated by reference in its entirety) are forms of amplification. Additional types of amplification include, but are not limited to, allele-specific PCR (see, e.g., U.S. Pat. No. 5,639,611; herein incorporated by reference in its entirety), assembly PCR (see, e.g., U.S. Pat. No. 5,965,408; herein incorporated by reference in its entirety), helicase-dependent amplification (see, e.g., U.S. Pat. No. 7,662,594; herein incorporated by reference in its entirety), Hot-start PCR (see, e.g., U.S. Pat. Nos. 5,773,258 and 5,338,671; each herein incorporated by reference in their entireties), intersequence-specfic PCR, inverse PCR (see, e.g., Triglia, et al et al. (1988) Nucleic Acids Res., 16:8186; herein incorporated by reference in its entirety), ligation-mediated PCR (see, e.g., Guilfoyle, R. et al et al., Nucleic Acids Research, 25:1854-1858 (1997); U.S. Pat. No. 5,508,169; each of which are herein incorporated by reference in their entireties), methylation-specific PCR (see, e.g., Herman, et al., (1996) PNAS 93(13) 9821-9826; herein incorporated by reference in its entirety), miniprimer PCR, multiplex ligation-dependent probe amplification (see, e.g., Schouten, et al., (2002) Nucleic Acids Research 30(12): e57; herein incorporated by reference in its entirety), multiplex PCR (see, e.g., Chamberlain, et al., (1988) Nucleic Acids Research 16(23) 11141-11156; Ballabio, et al., (1990) Human Genetics 84(6) 571-573; Hayden, et al., (2008) BMC Genetics 9:80; each of which are herein incorporated by reference in their entireties), nested PCR, overlap-extension PCR (see, e.g., Higuchi, et al., (1988) Nucleic Acids Research 16(15) 7351-7367; herein incorporated by reference in its entirety), real time PCR (see, e.g., Higuchi, et al et al., (1992) Biotechnology 10:413-417; Higuchi, et al., (1993) Biotechnology 11:1026-1030; each of which are herein incorporated by reference in their entireties), reverse transcription PCR (see, e.g., Bustin, S. A. (2000) J. Molecular Endocrinology 25:169-193; herein incorporated by reference in its entirety), solid phase PCR, thermal asymmetric interlaced PCR, and Touchdown PCR (see, e.g., Don, et al., Nucleic Acids Research (1991) 19(14) 4008; Roux, K. (1994) Biotechniques 16(5) 812-814; Hecker, et al., (1996) Biotechniques 20(3) 478-485; each of which are herein incorporated by reference in their entireties). Polynucleotide amplification also can be accomplished using digital PCR (see, e.g., Kalinina, et al., Nucleic Acids Research. 25; 1999-2004, (1997); Vogelstein and Kinzler, Proc Natl Acad Sci USA. 96; 9236-41, (1999); International Patent Publication No. WO05023091A2; US Patent Application Publication No. 20070202525; each of which are incorporated herein by reference in their entireties).

The term “polymerase chain reaction” (“PCR”) refers to the method of K. B. Mullis U.S. Pat. Nos. 4,683,195, 4,683,202, and 4,965,188, that describe a method for increasing the concentration of a segment of a target sequence in a mixture of genomic DNA without cloning or purification. This process for amplifying the target sequence consists of introducing a large excess of two oligonucleotide primers to the DNA mixture containing the desired target sequence, followed by a precise sequence of thermal cycling in the presence of a DNA polymerase. The two primers are complementary to their respective strands of the double stranded target sequence. To effect amplification, the mixture is denatured and the primers then annealed to their complementary sequences within the target molecule. Following annealing, the primers are extended with a polymerase so as to form a new pair of complementary strands. The steps of denaturation, primer annealing, and polymerase extension can be repeated many times (i.e., denaturation, annealing and extension constitute one “cycle”; there can be numerous “cycles”) to obtain a high concentration of an amplified segment of the desired target sequence. The length of the amplified segment of the desired target sequence is determined by the relative positions of the primers with respect to each other, and therefore, this length is a controllable parameter. By virtue of the repeating aspect of the process, the method is referred to as the “polymerase chain reaction” (“PCR”). Because the desired amplified segments of the target sequence become the predominant sequences (in terms of concentration) in the mixture, they are said to be “PCR amplified” and are “PCR products” or “amplicons.”

As used herein, the term “nucleic acid detection assay” refers to any method of determining the nucleotide composition of a nucleic acid of interest. Nucleic acid detection assay include but are not limited to, DNA sequencing methods, probe hybridization methods, structure specific cleavage assays (e.g., the INVADER assay, (Hologic, Inc.) and are described, e.g., in U.S. Pat. Nos. 5,846,717, 5,985,557, 5,994,069, 6,001,567, 6,090,543, and 6,872,816; Lyamichev et al., Nat. Biotech., 17:292 (1999), Hall et al., PNAS, USA, 97:8272 (2000), and US 2009/0253142, each of which is herein incorporated by reference in its entirety for all purposes); enzyme mismatch cleavage methods (e.g., Variagenics, U.S. Pat. Nos. 6,110,684, 5,958,692, 5,851,770, herein incorporated by reference in their entireties); polymerase chain reaction; branched hybridization methods (e.g., Chiron, U.S. Pat. Nos. 5,849,481, 5,710,264, 5,124,246, and 5,624,802, herein incorporated by reference in their entireties); rolling circle replication (e.g., U.S. Pat. Nos. 6,210,884, 6,183,960 and 6,235,502, herein incorporated by reference in their entireties); NASBA (e.g., U.S. Pat. No. 5,409,818, herein incorporated by reference in its entirety); molecular beacon technology (e.g., U.S. Pat. No. 6,150,097, herein incorporated by reference in its entirety); E-sensor technology (Motorola, U.S. Pat. Nos. 6,248,229, 6,221,583, 6,013,170, and 6,063,573, herein incorporated by reference in their entireties); cycling probe technology (e.g., U.S. Pat. Nos. 5,403,711, 5,011,769, and 5,660,988, herein incorporated by reference in their entireties); Dade Behring signal amplification methods (e.g., U.S. Pat. Nos. 6,121,001, 6,110,677, 5,914,230, 5,882,867, and 5,792,614, herein incorporated by reference in their entireties); ligase chain reaction (e.g., Barnay Proc. Natl. Acad. Sci USA 88, 189-93 (1991)); and sandwich hybridization methods (e.g., U.S. Pat. No. 5,288,609, herein incorporated by reference in its entirety).

As used herein, the terms “complementary” or “complementarity” used in reference to polynucleotides (i.e., a sequence of nucleotides) refers to polynucleotides related by the base-pairing rules. For example, the sequence “5′-A-G-T-3′,” is complementary to the sequence “3′-T-C-A-5′.” Complementarity may be “partial,” in which only some of the nucleic acids' bases are matched according to the base pairing rules. Or, there may be “complete” or “total” complementarity between the nucleic acids. The degree of complementarity between nucleic acid strands has significant effects on the efficiency and strength of hybridization between nucleic acid strands. This is of particular importance in amplification reactions, as well as detection methods that depend upon binding between nucleic acids.

As used herein, the term “primer” refers to an oligonucleotide, whether occurring naturally, as in a purified restriction digest, or produced synthetically, that is capable of acting as a point of initiation of synthesis when placed under conditions in which synthesis of a primer extension product that is complementary to a nucleic acid strand is induced (e.g., in the presence of nucleotides and an inducing agent such as a biocatalyst (e.g., a DNA polymerase or the like). The primer is typically single stranded for maximum efficiency in amplification, but may alternatively be partially or completely double stranded. The portion of the primer that hybridizes to a template nucleic acid is sufficiently long to prime the synthesis of extension products in the presence of the inducing agent. The exact lengths of the primers will depend on many factors, including temperature, source of primer and the use of the method. Primers may comprise labels, tags, capture moieties, etc.

As used herein, the term “nucleic acid molecule” refers to any nucleic acid containing molecule, including but not limited to, DNA or RNA. The term encompasses sequences that include any of the known base analogs of DNA and RNA including, but not limited to, 4 acetylcytosine, 8-hydroxy-N6-methyladenosine, aziridinylcytosine, pseudoisocytosine, 5-(carboxyhydroxyl-methyl) uracil, 5-fluorouracil, 5-bromouracil, 5-carboxymethylaminomethyl-2-thiouracil, 5-carboxymethyl-aminomethyluracil, dihydrouracil, inosine, N6-isopentenyladenine, 1-methyladenine, 1-methylpseudo-uracil, 1-methylguanine, 1-methylinosine, 2,2-dimethyl-guanine, 2-methyladenine, 2-methylguanine, 3-methyl-cytosine, 5-methylcytosine, N6-methyladenine, 7-methylguanine, 5-methylaminomethyluracil, 5-methoxy-amino-methyl-2-thiouracil, beta-D-mannosylqueosine, 5′-methoxycarbonylmethyluracil, 5-methoxyuracil, 2-methylthio-N-isopentenyladenine, uracil-5-oxyacetic acid methylester, uracil-5-oxyacetic acid, oxybutoxosine, pseudouracil, queosine, 2-thiocytosine, 5-methyl-2-thiouracil, 2-thiouracil, 4-thiouracil, 5-methyluracil, N-uracil-5-oxyacetic acid methylester, uracil-5-oxyacetic acid, pseudouracil, queosine, 2-thiocytosine, and 2,6-diaminopurine.

As used herein, the term “nucleobase” is synonymous with other terms in use in the art including “nucleotide,” “deoxynucleotide,” “nucleotide residue,” “deoxynucleotide residue,” “nucleotide triphosphate (NTP),” or deoxynucleotide triphosphate (dNTP).

An “oligonucleotide” refers to a nucleic acid that includes at least two nucleic acid monomer units (e.g., nucleotides), typically more than three monomer units, and more typically greater than ten monomer units. The exact size of an oligonucleotide generally depends on various factors, including the ultimate function or use of the oligonucleotide. To further illustrate, oligonucleotides are typically less than 200 residues long (e.g., between 15 and 100), however, as used herein, the term is also intended to encompass longer polynucleotide chains. Oligonucleotides are often referred to by their length. For example a 24 residue oligonucleotide is referred to as a “24-mer”. Typically, the nucleoside monomers are linked by phosphodiester bonds or analogs thereof, including phosphorothioate, phosphorodithioate, phosphoroselenoate, phosphorodiselenoate, phosphoroanilothioate, phosphoranilidate, phosphoramidate, and the like, including associated counterions, e.g., H⁺, NH₄ ⁺, Na⁺, and the like, if such counterions are present. Further, oligonucleotides are typically single-stranded. Oligonucleotides are optionally prepared by any suitable method, including, but not limited to, isolation of an existing or natural sequence, DNA replication or amplification, reverse transcription, cloning and restriction digestion of appropriate sequences, or direct chemical synthesis by a method such as the phosphotriester method of Narang et al. (1979) Meth Enzymol. 68: 90-99; the phosphodiester method of Brown et al. (1979) Meth Enzymol. 68: 109-151; the diethylphosphoramidite method of Beaucage et al. (1981) Tetrahedron Lett. 22: 1859-1862; the triester method of Matteucci et al. (1981) J Am Chem Soc. 103:3185-3191; automated synthesis methods; or the solid support method of U.S. Pat. No. 4,458,066, entitled “PROCESS FOR PREPARING POLYNUCLEOTIDES,” issued Jul. 3, 1984 to Caruthers et al., or other methods known to those skilled in the art. All of these references are incorporated by reference.

A “sequence” of a biopolymer refers to the order and identity of monomer units (e.g., nucleotides, amino acids, etc.) in the biopolymer. The sequence (e.g., base sequence) of a nucleic acid is typically read in the 5′ to 3′ direction.

The term “wild-type” refers to a gene or gene product that has the characteristics of that gene or gene product when isolated from a naturally occurring source. A wild-type gene is that which is most frequently observed in a population and is thus arbitrarily designed the “normal” or “wild-type” form of the gene. In contrast, the terms “modified,” “mutant,” and “variant” refer to a gene or gene product that displays modifications in sequence and or functional properties (i.e., altered characteristics) when compared to the wild-type gene or gene product. It is noted that naturally occurring mutants can be isolated; these are identified by the fact that they have altered characteristics when compared to the wild-type gene or gene product.

As used herein, the term “gene” refers to a nucleic acid (e.g., DNA) sequence that comprises coding sequences necessary for the production of a polypeptide, precursor, or RNA (e.g., rRNA, tRNA). The polypeptide can be encoded by a full length coding sequence or by any portion of the coding sequence so long as the desired activity or functional properties (e.g., enzymatic activity, ligand binding, signal transduction, immunogenicity, etc.) of the full-length or fragment polypeptide are retained. The term also encompasses the coding region of a structural gene and the sequences located adjacent to the coding region on both the 5′ and 3′ ends for a distance of about 1 kb or more on either end such that the gene corresponds to the length of the full-length mRNA. Sequences located 5′ of the coding region and present on the mRNA are referred to as 5′ non-translated sequences. Sequences located 3′ or downstream of the coding region and present on the mRNA are referred to as 3′ non-translated sequences. The term “gene” encompasses both cDNA and genomic forms of a gene. A genomic form or clone of a gene contains the coding region interrupted with non-coding sequences termed “introns” or “intervening regions” or “intervening sequences.” Introns are segments of a gene that are transcribed into nuclear RNA (e.g., hnRNA); introns may contain regulatory elements (e.g., enhancers). Introns are removed or “spliced out” from the nuclear or primary transcript; introns therefore are absent in the messenger RNA (mRNA) transcript. The mRNA functions during translation to specify the sequence or order of amino acids in a nascent polypeptide.

In addition to containing introns, genomic forms of a gene may also include sequences located on both the 5′ and 3′ end of the sequences that are present on the RNA transcript. These sequences are referred to as “flanking” sequences or regions (these flanking sequences are located 5′ or 3′ to the non-translated sequences present on the mRNA transcript). The 5′ flanking region may contain regulatory sequences such as promoters and enhancers that control or influence the transcription of the gene. The 3′ flanking region may contain sequences that direct the termination of transcription, post-transcriptional cleavage and polyadenylation.

As used herein, the terms “multimethylation,” “series methylation” and “specific methylation” are used interchangeably to refer to defined combinations of CpG sites or loci in a marker sequence must be methylated to call that sequence methylated in a coordinate or multimethylation assay. For example, a specific methylation assay of the CpG sites for BMP3 might require that the CpG positions at 23, 34, 53, 61, 70, and 74, numbered by reference to FIGS. 1A and 1B, all be methylated in order for a sample to be classified as methylated at the BMP3 marker. Specific methylation of BMP3 is not limited to this set of particular loci, but may include more, fewer, or a different collection of CpG loci. The CpG loci selected for co-analysis in a multimethylation assay are preferably selected. e.g., by analysis of normal (non-adenoma, non-cancer) samples to identify CpG methylation combinations that are less frequently represented in normal samples. In preferred embodiments, combinations of methylation sites are selected to produce good signal-to-noise in cancer and adenoma samples (i.e., the mean multimethylation at a particular combination of loci in cancer samples divided by the mean multimethylation in at those loci in normal samples is high).

As used herein, the terms “individual” and “average” methylation are used interchangeably to refer to analyses in which each CpG locus is analyzed individually, such that all molecules in which that base is methylated are included in a count, regardless of the methylation status of other loci, e.g., in the same marker. Generally, the methylation percentages of all the loci in a marker/region are then averaged, to produce a percent methylation figure for that marker.

As used herein, the term “kit” refers to any delivery system for delivering materials. In the context of reaction assays, such delivery systems include systems that allow for the storage, transport, or delivery of reaction reagents (e.g., oligonucleotides, enzymes, etc. in the appropriate containers) and/or supporting materials (e.g., buffers, written instructions for performing the assay etc.) from one location to another. For example, kits include one or more enclosures (e.g., boxes) containing the relevant reaction reagents and/or supporting materials. As used herein, the term “fragmented kit” refers to a delivery systems comprising two or more separate containers that each contain a subportion of the total kit components. The containers may be delivered to the intended recipient together or separately. For example, a first container may contain an enzyme for use in an assay, while a second container contains oligonucleotides. The term “fragmented kit” is intended to encompass kits containing Analyte specific reagents (ASR's) regulated under section 520(e) of the Federal Food, Drug, and Cosmetic Act, but are not limited thereto. Indeed, any delivery system comprising two or more separate containers that each contains a subportion of the total kit components are included in the term “fragmented kit.” In contrast, a “combined kit” refers to a delivery system containing all of the components of a reaction assay in a single container (e.g., in a single box housing each of the desired components). The term “kit” includes both fragmented and combined kits.

As used herein, the term “information” refers to any collection of facts or data. In reference to information stored or processed using a computer system(s), including but not limited to internets, the term refers to any data stored in any format (e.g., analog, digital, optical, etc.). As used herein, the term “information related to a subject” refers to facts or data pertaining to a subject (e.g., a human, plant, or animal). The term “genomic information” refers to information pertaining to a genome including, but not limited to, nucleic acid sequences, genes, allele frequencies, RNA expression levels, protein expression, phenotypes correlating to genotypes, etc. “Allele frequency information” refers to facts or data pertaining to allele frequencies, including, but not limited to, allele identities, statistical correlations between the presence of an allele and a characteristic of a subject (e.g., a human subject), the presence or absence of an allele in an individual or population, the percentage likelihood of an allele being present in an individual having one or more particular characteristics, etc.

DESCRIPTION OF THE DRAWINGS

FIGS. 1A and 1B provide sequence and CpG information for exemplary marker regions used in the present analysis. For each target gene, the native sequence of the region is shown in the top line. Unmethylated C-residues that would be converted by bisulfite and amplification to Ts are shown as T residues. Candidate methylation positions are shown boxed. Reference numbering for base and CpG positions is shown above each native sequence. Primer locations for amplification are shown as a row of underlined base positions.

FIGS. 2 A-J provide tables showing analyses of normal, adenoma, and cancer samples in which the average methylation was determined at each of the indicated CpG positions, in the indicated marker regions. For each marker, the numbered CpG positions are as indicated by the reference numbers in FIGS. 1A and 1B. The Mean methylation at each specific locus is shown at the bottom of each column for normal, adenoma, and cancer samples. The ratio of normal/mutant methylation for each locus (a methylation ratio at each locus) is shown at the bottom of each column of the Adenoma and Cancer sample data. The Mean columns on the right of each table indicate the average of methylation across all indicated CpG loci for each of the samples. The Mean and SD values across all normal samples at all loci are as indicated below each table of values from normal samples.

FIGS. 3 A-I provide tables showing the analyses of normal, adenoma, and cancer samples in which average methylation across all of the CpG loci indicated in FIGS. 2A-J were calculated for each marker in each sample. For the normal samples in FIG. 3A, the average, standard deviation and the mean plus 2 or 3 standard deviations for each marker are indicated. For the adenoma and cancer samples, shaded cells in FIGS. 3B and 3C indicate a positive result, reflected as an average methylation value for that marker that is greater than the mean methylation+3 standard deviations determined for that marker in the normal samples.

FIGS. 3D and 3E show the calculated effect of a 20-fold dilution of adenoma and cancer DNA into normal DNA, FIGS. 3F and 3G show a calculated 10-fold dilution, and FIGS. 3H and 3I show a calculated 5-fold dilution. In each of the calculated dilutions, the average methylation for a marker is divided by the 20, 10, or 5, added to the mean methylation of the normal DNA for that marker. Shaded cells in FIGS. 3D-3I indicate an average methylation value for that marker that is greater than the mean methylation+2 standard deviations (specificity of 97.5%) determined for that marker in the normal samples.

Below each of FIGS. 3B-3I, the percentage of positive values for each marker in the sample type and dilution for that panel are indicated. The percentage of samples giving a positive signal for at least one of the Vimentin, BMP3, Septin 9 and TFPI2 markers are indicated at the bottom of each panel.

FIGS. 4A and 4B provide sequence and CpG information for exemplary genes used in the present analysis. The CpG loci in each marker gene included in the defined subsets of CpG loci for coordinate methylation analysis in colorectal adenoma and cancer samples are shown with a black background and in white typeface.

FIGS. 5 A-I provide tables showing the analyses of normal, adenoma, and cancer samples in which methylation was determined at each of the indicated CpG positions in the indicated marker regions (i.e., samples were assayed for the percentage of DNA copies that displayed methylation at all of the CpG loci in the defined subset). Each marker was tested at each of the CpG loci in the defined subsets indicated in FIGS. 4A and 4B and the percentage methylation data reflects the percentage of marker copies having methylation at all of the tested CpG loci (coordinate methylation or “multimethylation” analysis). For the normal samples in FIG. 5A, the mean multimethylation, standard deviation and the mean plus 2 or 3 standard deviations for each marker are indicated. For the adenoma and cancer samples, shaded cells in FIGS. 5B and 5C indicate a positive result, reflected as multimethylation value for that marker that is greater than the mean multimethylation+3 standard deviations determined for that marker in the normal samples.

FIGS. 5D and 5E show the calculated effect of a 20-fold dilution of adenoma and cancer DNA into normal DNA, FIGS. 5F and 5G show a calculated 10-fold dilution, and 5H and 5I show a calculated 5-fold dilution. In each of the calculated dilutions, the average multimethylation for a marker is divide by the 20, 10, or 5, added to the mean multimethylation of the normal DNA for that marker. Shaded cells in FIGS. 5D-5I indicate an average multimethylation value for that marker that is greater than the mean multimethylation+2 standard deviations (specificity of 97.5%) determined for that marker in the normal samples.

Below each of FIGS. 5B-5I, the percentage of positive values for each marker in the sample type and dilution for that panel are indicated. The percentage of samples giving a positive signal for at least one of the vimentin, BMP3, Septin 9 and TFPI2 markers are indicated at the bottom of each panel.

FIG. 6 shows a table and graph comparing the percent positive values calculated for each marker in adenoma and cancer samples, as indicated, using either individual/average methylation or multimethylation analysis methods to test each of the indicated markers, at each of the indicated calculated dilutions.

FIG. 7 shows a table and graph comparing the percent positive values determined in adenoma and cancer samples, determined using the four markers with the lowest mean background in these samples (vimentin, BMP3, Septin9, TFPI2), using either the individual/average methylation or the multimethylation analysis method, at each of the indicated calculated dilutions into normal DNA.

DETAILED DESCRIPTION OF THE INVENTION

Embodiments of the invention are described in this summary, and in the Summary of the Invention, above, which is incorporated here by reference. Although the invention has been described in connection with specific embodiments, it should be understood that the invention as claimed should not be unduly limited to such specific embodiments.

The present invention relates to methods and compositions for determination of, and uses of, specific methylation patterns indicative of adenoma and carcinoma. In particular, the invention relates to analysis of defined subsets of CpG loci that are coordinately methylated in DNAs from cancer and adenoma samples, methods for identifying coordinately methylated loci, and methods of using analysis of coordinately methylated loci in one or more markers or regions in the design of assays for adenoma and cancer having improved sensitivity and specificity.

The present invention relates to the observation that, within marker nucleic acids for which methylation status is indicative of cellular status, e.g., cancerous, pre-cancerous, normal, etc., a subset of the individual methylation loci, e.g., CpG loci, in non-normal cells generally displays a greater degree of methylation relative to the background levels of methylation observed at corresponding loci in normal cells, while other methylation loci in the non-normal cells may exhibit levels of methylation that are closer to background levels. In some embodiments, the degree of methylation observed for a particular locus in plurality of cancerous or pre-cancerous cells relative to normal cells is expressed as a methylation ratio.

Some embodiments of the present invention relate to screening known or suspected marker genes to identify specific methylation loci that exhibit greater ratios of disease-associated methylation relative to background methylation, as compared to other marker genes or other loci in the same marker gene. In some preferred embodiments, the present invention relates to coordinate methylation analysis, to measure the degree to which a marker molecule or a sample exhibits methylation at all of a plurality of selected loci.

The present invention relates to analyzing methylation statuses of a defined set of individual CpG loci in methylation markers (or target regions within such markers) in a significant enough number of individual DNA molecules in adenoma samples or cancer samples to identify defined subsets of CpG loci that have advantageous methylation ratios compared to other loci in the same adenoma or cancer samples. A defined subset of CpG loci that have advantageous methylation ratios in a sample may comprise the entirety of a set of CpG loci in a particular marker or target region of a marker, or it may be fewer than all of the CpG loci in the characterized region of the marker.

Conventional methods of analyzing methylation status of a marker generally involve analysis of a mixed population of molecules. For example, amplification of a marker nucleic acid from a sample generally produces a mixture of amplicons coming from many copies of a target molecule. If the amplification conditions are not selective for a gene variant, the amplicon product contains a mixture of the variant and the normal or wildtype DNA. Even if primers are specific for a mutation or for a particular methylation site, when DNA is amplified from many copies of target DNA derived from many cells, there can be heterogeneity in other base positions in the resulting amplicon. If these mixed amplicons are sequenced directly, the resulting sequence traces reveal the consensus sequence of the mixed population, and particular sequences or mutations present in a small portion of the population are essentially undetectable. Although some researchers have sequenced individual clones from such amplifications to examine sequence information for individual molecules from the mixture, only small numbers of molecules have been analyzed and the data gathered did not suggest any that any specific loci within the markers predictably exhibited advantageous methylation ratios compared to other methylated loci within the same target, or that coordindated analysis of loci having advantageous methylation ratios could be useful in improving the specificity and sensitivity of assays in detection of neoplasms. An aspect of the present invention is based on the observation that collecting methylation ratio information from a very large number of individual molecules in both normal and non-normal samples reveals that some methylation loci in marker regions or sequences exhibit a greater degree of methylation in non-normal cells compared to background (methylation at the same loci in normal cells) than do other individual loci in the same marker region or gene. These loci in non-normal sequences that have a greater level of methylation compared to background can be viewed as being particularly advantageous in that they are easier to identify over the background level of methylation observed in normal cells. One aspect of this advantage is that analysis of these particular loci permits identification of cancer-associated methylation with more sensitivity, and in a greater background of normal cells.

The present invention also relates to the observation that coordinated analysis of multiple loci provides a significantly enhanced level of sensitivity in the identification of cancerous or precancerous cells, especially in samples that may also comprise a significant number of normal cells. For example, FIG. 6 compares the sensitivity of detecting adenoma and cancer cells. For each of the indicated marker genes, the methylation was either determined as an average across the marker region (e.g., the mean methylation in the vimentin marker across all of loci 26, 37, 40, 45, 52, 54, 59, 63, and 74; see FIGS. 3A-I), indicated as “Individual” average methylation, or as a percentage of molecules displaying methylation at all of a subset of selected loci (e.g., methylation in the vimentin marker at all three of loci 37, 40 and 45; see FIGS. 5A-I), i.e., coordinate methylation analysis of multiple individual loci, indicated as “multi”. The sensitivities for the same samples are also shown in calculated 5, 10 or 20-fold dilutions into normal DNA. FIG. 6 shows that, while assay sensitivities may be similar in DNA analyzed directly from tissue without dilution, as the amount background from normal DNA is increased at the larger dilutions, the coordinate methylation analysis is shown to be far more sensitive than the average methylation analysis. For example, in the analysis of Septin9, the adenoma and cancer samples can be detected above background only in the undiluted and 5-fold dilution profiles when average methylation across the marker is analyzed, while these same samples can be detected with over about 69-74% sensitivity at 20-fold dilution, and 90-93% sensitivity at 10-fold dilution, when coordinate methylation analysis of loci 37, 40 and 45 is used.

In some embodiments, the present invention provides a method for designing a methylation assay to identify a disease state, comprising I) selecting at least one sequence for analysis; II) determining the methylation status of a plurality of loci in the sequence in a population of normal cells and a population of non-normal cells to determine an average rate of methylation for each of the plurality of loci each both normal and non-normal cells; and III) identifying at least two loci in said plurality of loci having advantageous methylation ratios.

I. Selection of Sequence(s).

Methylation markers associated with particular disease states have been identified for a number of disease states. For example, colorectal neoplasm-specific marker include, e.g., bmp-3, bmp-4, SFRP2, vimentin, septin9, ALX4, EYA4, TFPI2, NDRG4, FOXE1, long DNA, BAT-26, K-ras, APC, melanoma antigen gene, p53, BRAF, and PIK3CA. Additional markers include but are not limited those in Table 1, above. Analysis of candidate methylation loci to identify those with advantageous methylation ratios may comprise analysis of every locus in a target sequence (e.g., every CpG) or it may comprise analysis of a subset of the methylation loci. In some embodiments, CpGs are selected for analysis by their location in particular methylation hotspots, while in other embodiments, loci for analysis may be conveniently located with respect to primer binding sites or other sequence features. FIG. 1A provides an exemplary selection of marker neoplasm-associated markers with each of the C residues in the CpG loci indicated by a box. For each target gene, the native sequence of the region is shown in the top line and the sequences for unmethylated and methulated DNA as they would appear following bisulfite conversion and amplification are shown below. Unmethylated C-residues that would be converted by bisulfite and amplification to T residues are shown as Ts.

In some embodiments, the present invention provides use of a nucleic acid detection assay to coordinately analyze a plurality of the advantageous loci in a sample, thereby determining the disease state of cells in sample.

II. Determining Methylation Ratios for Loci in the Selected Sequence(s).

As discussed above, determining a methylation ratio for a locus comprises determining an average rate of methylation for that locus in a population of normal cells and determining the average rate of methylation at the same locus in a population of non-normal cells. As noted above, commonly used methods of methylation analysis of marker genes are performed on mixed nucleic acids, e.g., amplicons produced from unfractionated DNA from a mixed population of cells (such as DNA purified from a multicellular tissue sample). While some studies have analyzed individual clones of amplicons made from unfractionated sample DNA, the numbers of clones analyzed has typically been too small to reveal significant or reproducible differences in methylation ratios at individual CpG loci within the sequences. For example, in their comparison of highly rmethylated genes in colorectal cancer, Zou, et al., analyzed only six clones from each sample. (Zou, et al., Cancer Epidemiol Biomarkers Prev 2007; 16(12):2686), while Weisenberg, et al. used The present invention comprises large-scale analysis of individual DNA molecules, for example, by direct sequencing of individual DNA molecules, or by sequencing of clonally amplified DNA.

While not limiting the present invention to any particular methods, methods of clonally amplifying individual copies of nucleic acids (e.g., using PCR) can be used in the rapid analysis of large numbers of individual markers from normal and non-normal samples. Single-molecule amplification methods may comprise use of microchambers, emulsion reactions, “bridge PCR” on solid supports, or any of a number of established methods of segregating the amplification products arising from individual target molecules. Following single molecule amplification, amplicons can be sequenced.

Improved methods of sequencing individual molecules directly obviate the need to clone molecules into cells or, in some methods, the need to clonally amplify prior to sequencing. Elimination of cloning into cells makes analysis of much larger collections of molecules dramatically more efficient. Platforms for individual molecule sequencing include the 454 FLX™ or 454 TITANIUM™ (Roche), the SOLEXA™/Illumina Genome Analyzer (Illumina), the HELISCOPE™ Single Molecule Sequencer (Helicos Biosciences), the Ion Personal Genome Machine (Ion Torrent), and the SOLID™ DNA Sequencer (Life Technologies/Applied Biosystems) instruments, as well as other platforms still under development by companies such as Intelligent Biosystems and Pacific Biosystems. Although the chemistry by which sequence information is generated varies for the different next-generation sequencing platforms, all of them share the common feature of generating sequence data from a very large number of individual sequencing templates, in sequencing reactions that are run simultaneously. Data from the reactions are collected using, e.g., a flow cell, a chemical or optical sensor, and/or scanner, and sequences are assembled and analyzed using bioinformatics software.

In certain preferred embodiments, the present invention provides methods of analysis of methylation markers using digital sequencing to identify neoplasm-associated methylation loci that have methylation ratios that are statistically significantly advantageous compared to other loci in the same markers. In preferred embodiments, digital sequencing is done in a highly or massively parallel fashion, providing higher precision in identifying CpG methylation sites having advantageous methylation ratios.

For the massively parallel digital sequencing methods mentioned above, each molecule is analyzed for methylation at each CpG locus, so the percentage of DNA copies having methylation at any combination of the CpG loci can be analyzed after the experimental run. Further, each particular marker sequence, e.g., each target nucleic acid molecule, or clonal amplicon may be interrogated many, many times, e.g., at least 100 times, sometimes over 1000 times, and in some instances over 100,000 times, or as many as 500,000 times. Thus, patterns of coordinate methylation indicative of cancer or adenoma that would be undetectable in analysis of a handful of individual target molecules may be revealed.

III. Selection of a Subset of Methylation Loci for Coordinate Analysis

As noted above, determination of the methylation status of a set of CpG loci in a large number of copies marker DNA from both normal samples and non-normal samples (e.g., adenoma or cancer samples) reveals that certain certain CpG loci in marker genes or regions may tend to be coordinately methylated. Further, design of nucleic acid detection assays to interrogate a plurality of CpG loci for which coordinate methylation is indicative of adenoma or cancer in a sample can provide an assay that has improved signal-to-noise compared to assays that survey average percent methylation across entire marker genes.

One aspect of selecting a subset of CpG loci comprises selecting loci that have been determined to be coordinately methylated by use, e.g., of digital analysis methods. Another aspect comprises selecting CpG loci determined to have advantageous methylation ratios when normal DNA is compared to adenoma or cancer DNA. Assay designs may, but need not, make use of a CpG locus having the most advantageous methylation ratio compared to other loci in the same marker. In some embodiments, selection of a plurality of CpG loci as a subset comprises selecting the plurality of loci having the most advantageous methylation ratios. In other embodiments, selection of a plurality of CpG loci as a subset comprises selecting the locus having the most advantageous methylation ratio, then selecting at least additional CpG loci that are conveniently situated with respect to the first selected locus for the configuration of a particular nucleic acid detection assay (e.g., the selection of CpG loci having particular proximity to each other for configuring an invasive cleavage assay, ligation assay, amplification assay, etc.) in order to interrogate all of the selected CpG loci on copies of the target DNA in a single assay. In some embodiments, a candidate subset of CpG loci is further analyzed to determine the percentage of copies of marker DNA from non-normal samples that are coordinately methylated at those candidate loci, and that have little or no coordinated methylation in normal samples.

Analysis of Samples for Detection of Adenoma or Cancer

Conventional methods of methylation analysis, (e.g., conventional methylation-specific PCR, real time methylation-specific PCR, see, e.g., U.S. Pat. Nos. 5,786,146, 6,017,704, 6,200,756, 6,265,171,), typically analyze in a non-digital fashion, e.g., analyzing a mixture of co-amplified molecules derived from a mixture of DNA target nucleic acids, so that analysis of the amplified products provides sequence information that reflects that aggregate or average methylation status in the amplicon population, but does not provide information on the percentage of starting molecules having coordinated methylation at all of a plurality of CpG loci. In some instances researchers have analyzed a number of cloned amplicons, which can reveal the diversity in methylation in the CpG loci within a target marker gene. However, sequencing individual clones has not provided enough data to reveal statistically significant coordinate methylation of specific subsets of CpG loci.

In contrast to conventional methods, we sought to analyze methylation marker genes in a massively parallel digital sequencing fashion to identify statistically significant coordinate methylation of specific CpG loci associated with neoplasms (adenoma and carcinoma). This method of analysis allowed us to:

1. Analyze samples for coordinate methylation in a marker gene as a means of detecting neoplasms without the need for testing any genetic (mutation) markers 2. Analyze samples for coordinate methylation in a plurality of marker genes as a means of detecting neoplasms without the need for testing any genetic (mutation) markers

We decided to use “digital” sequencing on a larger number of tissue samples obtained by biopsy from colorectal adenomas, colorectal cancers normal colorectal epithelia and other GI cancers and sequence a number of specific regions within several genes. This type of sequencing provides a methylation pattern for each individual methylated gene. For the first run we have 9 normal tissues, 38 adenomas and 36 cancer samples with the following markers—Vimentin, BMP3, Septin 9, TFPI2, 2 regions of LRAT, and EYA4.

Surprisingly we found that in some of the genes, the background observed as methylation in normal samples is randomly distributed in the sequences, while the methylation associated with cancer and adenoma is not. Thus, if certain rules are applied e.g. if all of C residues a, b, and c have to methylated in a diagnostic assay, then the number of DNA copies presenting methylation at all three positions in that sequence is reduced compared to the number of DNA copies displaying methylation at a subset of the positions. In some of the marker genes or regions tested, the reduction in number of DNA copies in normal DNA displaying methylation at all of the selected sites drops to a greater degree than it does in the DNA from cancer and/or adenoma sample, resulting in an significantly enhanced ratio of specific signal to background noise. For certain genes, the background from normal DNA is dramatically reduced by using multimethylation (coordinate methylation) analysis, while no equivalent reduction in signal from cancer and adenoma DNA is seen. For other genes the background in normal samples is less reduced and/or the signal from cancer DNA also decreases with multimethylation analysis, such that there is less or no net improvement in the signal-to-noise and the advantages of using a multimethylation analysis approached are less. Genes having favorable signal to noise in multimethylation analyses are readily determined empirically.

From the multimethyation data presented herein (see, e.g., FIGS. 2A-J) it is possible to:

-   -   a. Identify regions within the gene sequences that give higher         discrimination between normal and non-normal (e.g., cancer and         adenoma) cells;     -   b. Identify particular genes that have greater signal-to-noise         (non-normal cell signal compared to normal cell background);     -   c. Identify particular methylation loci that have greater         signal-to-noise;     -   d. Identify particular methylation loci that are coordinately         methylation in adenoma and cancer samples but not in normal         samples, such that detection coordinate methylation at these         loci is a sensitive indicator of adenoma or cancer;     -   e. Identify genes with very low background methylation allowing         for greater dilution of methylated DNA in normal DNA with less         decrease in assay sensitivity;     -   f. Identify genes that are diagnostically complementary to each         other, i.e., that when analyzed in combination produce         diagnostic information of elevated sensitivity and/or elevated         specificity compared to the genes analyzed alone.     -   g. Identify combinations of genes that give elevated sensitivity         at elevated specificity, e.g., 100% sensitivity for cancer and         adenoma at 100% specificity.

EXPERIMENTAL EXAMPLES Example 1 Use of Digital PCR and Sequencing for the Identification of Specific Subsets of CpG Loci Methylated in Cancer and Adenoma Samples

DNA extracted from frozen tissue samples was treated with an EPITECT bisulfite conversion kit (Qiagen) to convert non-methylated cytosines to uracil. Methylated cytosines remain unconverted. Primers for each gene region were designed for each sequence such that the composition of the amplification products remained the same as the original target sequences and methylated and non methylated sequences were amplified with equal efficiencies. Amplification of the dU-containing converted DNA produced amplicons having T-residues in place of the dU residues. The amplicons were then prepared for sequencing on the Illumina instrument. For each tissue sample, the amplification reaction for each target was prepared from the same sample of bisulfite-treated DNA.

After sequencing, the data was analyzed quantitatively as an average methylation similar to Sanger sequencing, but at a higher precision and resolution in that the combined signal at each position is calculated from individual molecules. For each amplicon sequence, a set of CpG loci was evaluated for percent methylation in the different tissues, to identify a subset the loci that were co-methylated more frequently in cancer and/or adenoma samples than in normal tissues.

Illumina Sequencing Protocol:

Sequencing was conducted according to the procedure recommended for the Illumina Genome Analyzer IIx, GAIIx, Data collection software ver. 2.5, and Pipeline analysis software ver 1.5. Briefly, the Illumina procedure comprises a) preparation of a library from sample DNA by attachment of known sequence tags that permit indexing, flow cell attachment, amplification, and sequencing; b) attachment of the library to a flow cell surface; c) bridge amplification to produce clusters of DNA fragments derived from single molecules, and d) sequencing in using iterative primer extension reactions using labeled reversible terminators to determine the nucleotide sequence of each cluster of amplicons. See, e.g., Bentley, et al, Nature 456, 53-59 (6 Nov. 2008)/doi:10.1038/nature07517 with supplementary methods and data, incorporated herein by reference. Use of unique tag sequences for indexing permits analysis of multiple samples in a single flow cell. See, e.g., Craig, et al., Nat. Methods Nat Methods. 2008 October; 5(10):887-93 (Epub 14 Sep. 2008), incorporated herein by reference.

Sample Set:

N=82, composed of tissue DNA extracted from 42 colorectal cancers, 31 pre-cancerous adenomas and 9 normal colonic mucosa.

Flow Cell Configuration

Samples were indexed at 12 per lane for a total of 7 lanes. A flow cell is composed of 8 lanes, one of which is dedicated to a phiX quality control.

Library Preparation:

Tissue-extracted DNA from patients was bisulfite-treated and a 2-step amplification using approximately 10,000 genome copies of initial material was carried out. The first round used tailed (T1)(Illumina) primers specific for marker sequences. These tails were Illumina-derived sequences needed for round two. The second round (T2)(Illumina) PCR uses primers specific for the Illmuna tails added in T1, and incorporates the index, sequencing primer, and flow cell attachment sequences. During the library preparation, multiple qPCR checks were run on the samples to ensure equimolar representations of all amplicons in the libraries.

Primer Design:

Forward and reverse primers specific for regions with converted, non-CpG cytosines are designed (using, e.g., MethPrimer software) to amplify each of the specific biomarker sites in a non-methylation specific manner. When CpG cytosines cannot be avoided in the primer design, degenerate mixtures (C/T; G/A) are used those sites in the primers. If additional sequences outside the primary amplicon need to be queried, additional primers may be designed. If CpGs in the target sequence cannot be avoided, the primers may incorporate degenerate bases at CpG sites (BiSearch software).

Primers for second round PCR comprise sequences for Illumina flow cell attachment (bridge amplification sites), sequencing primer sites (for the sample read), index sites, and sequencing primers sites (for the indexing read). Each of the primer sets (x) has 12 different index tags, for a total of 12× sets. Index-independent primer sets (n=x) are optimized on converted non-methylated DNA (e.g., human DNA) and converted methylated DNA. For example, the control DNAs are amplified, purified (e.g, using AMPURE treatment (Agencourt)), and run on an Agilent 2100 Bioanalyzer to assess the size and quantity of the amplified nucleic acids.

Experimental Steps: DNA Isolation and Bisulfite Conversion:

-   1) DNA was extracted and purified from tissue using either DNAZOL     (Invitrogen), or QIAAMP kits (Qiagen) and the concentration and     purity was measured by absorbance (A230/A260/A280) using a Nanodrop     ND-1000 spectrophotometer (Thermo Scientific). PICOGREEN     fluorescence (Molecular Probes) was used in conjunction with a TECAN     F-200 (Tecan) plate-reader for high samples exhibiting high A230     values. -   2) Samples were adjusted to a concentration of at least 200 ng/uL     using a Speedvac evaporation concentrator (Thermo Scientific), as     necessary. -   3) For each sample, 2 ug of DNA was bisulfite-treated using EPITECT     96-well plates (Qiagen). -   4) Recovery was assessed by absorbance and OLIGOGREEN fluorescence     (Molecular Probes) and the conversion efficiency was assessed with     quantitative PCR using cytosine-containing, non-CpG primers specific     for unconverted DNA. Conversion efficiency was determined to be     greater than 99%

First-Round PCR:

-   5) The 84 samples were amplified in reactions with 30 ng DNA using     marker-specific (T1) primer sets. The number of cycles used was     specific for each marker set and was empirically determined by the     initial control reactions on both methylated and unmethylated DNA.     The number of cycles is approximately set at the mean calculated Ct     value. For example, the following numbers of cycles were used for     the indicated markers:     -   TFPI12; 26 cycles     -   SEPT9; 27 cycles     -   BMP3; 28 cycles     -   VIM; 28 cycles     -   EYA4; 29 cycles -   6) The amplified product from each reaction was purified using     AMPURE beads (Agencourt), with elution in EB buffer (Qiagen). -   7) The product for each marker was quantified by qPCR as described     above, using a T2 primer set. Master plates were prepared containing     equal amounts of each biomarker for each sample.

Second-Round PCR:

-   8) The first-round samples were then amplified with the 12 T2     indexed primers. -   9) The product for each reaction was again purified and     concentrations measured with qPCR, this time with flow-cell-specific     primers and a standard curve created using serial dilution of PhiX     control DNA.

Final Library Prep:

-   10) The 12 columns of each plate were combined into 1 master column     in equimolar proportions. 1 uL of each library was loaded on a high     sensitivity DNA chip (Agilent) and run on the Bioanalyzer. A final     qPCR was also performed on the 480 LightCycler (Roche) with the PhiX     standards. -   11) The libraries were sequenced on the Illumina Instrument and the     sequence data obtained for each marker for each sample. -   12) Average methylation at each CpG site was calculated for each     marker for each sample. See FIGS. 2A-J. -   13) The percentage of molecules methylated at all of the CpG loci in     a defined subset of the CpG loci in each marker was calculated for     each marker for each sample. See FIGS. 5A-I.

All publications and patents mentioned in the above specification are herein incorporated by reference. Various modifications and variations of the described methods and systems of the invention will be apparent to those skilled in the art without departing from the scope and spirit of the invention. Although the invention has been described in connection with specific preferred embodiments, it should be understood that the invention as claimed should not be unduly limited to such specific embodiments. Indeed, various modifications of the described modes for carrying out the invention that are obvious to those skilled in the relevant fields are intended to be within the scope of the following claims. 

1-30. (canceled)
 31. A method of characterizing a sample comprising: i) providing a selected defined subset of CpG loci for at least one marker DNA in said sample, wherein said defined subset of CpG loci has been selected by a process comprising: a) isolating marker DNA from a first plurality of samples in which said marker DNA is unmethylated DNA that exhibits background methylation, and from a second plurality of samples in which said marker DNA is methylated; b) treating said isolated marker DNA with bisulfite under conditions wherein unmethylated cytosine residues are converted to uracil residues; c) determining the methylation status for three or more CpG loci in each of at least 1000 individual copies of said isolated marker DNA from each of said first plurality and said second plurality of samples; d) for each of said three or more CpG loci, determining a ratio between a mean methylation at that CpG locus in said first plurality of samples to a mean methylation at the corresponding CpG locus in said second plurality of samples; and e) selecting a defined subset of at least three CpG loci for which the percentage of individual copies of the marker DNA from said first plurality of samples that are methylated at all of said at least three CpG loci in said defined subset is less than the percentage of individual copies of said marker DNA from said second plurality of samples that are methylated at all of said at least three CpG loci in said defined subset; ii) providing bisulfite-treated marker DNA from said sample, and iii) characterizing the methylation state of said at least one marker DNA, wherein characterizing said the methylation state of said at least one marker DNA comprises detecting the methylation status of all of the at least three CpG loci in said defined subset of CpG loci in a strand of the bisulfite-treated marker DNA, wherein the absence of methylation at any of the CpG loci in said defined subset of CpG loci in a strand of said marker DNA produces an assay result classifying that strand of said marker DNA as not methylated; iv) reporting the assay result for said at least one marker DNA from said sample.
 32. The method of claim 31, wherein said defined subset of CpG loci comprises at least four CpG loci.
 33. The method of claim 31, wherein said defined subset of CpG loci comprises at least five CpG loci.
 34. The method of claim 31, wherein said determining comprises analysis of said CpG loci in a nucleic acid detection assay configured to determine the methylation status of each of said loci in a single nucleic acid detection assay.
 35. The method of claim 34, wherein said determining comprises analysis of said CpG loci in a nucleic acid detection assay configured to determine the methylation status of each of said loci in a single reaction mixture.
 36. The method of claim 31, wherein said DNA detection assay comprises at least one assay selected from the group consisting of a primer extension assay, a nucleic acid amplification assay, a structure-specific cleavage assay, 5′ nuclease cleavage assay, an invasive cleavage assay, and a ligation assay.
 37. The method of claim 31, wherein said sample is a sample from a human.
 38. The method of claim 37, wherein said sample from a human is a stool sample.
 39. The method of claim 37, wherein said at least one marker DNA comprises cancer and/or adenoma marker DNA.
 40. The method of claim 39, wherein said at least one marker DNA comprises at least three cancer or adenoma marker DNAs.
 41. The method of claim 39, wherein said at least one marker DNA comprises at least four cancer or adenoma marker DNAs.
 42. The method of claim 39, wherein said cancer or adenoma marker DNA is selected from the group consisting of vimentin, BMP3, Septin 9, TFPI2, 2 regions of LRAT, and EYA4 DNAs.
 43. The method of claim 42, wherein said at least one marker DNA comprises vimentin DNA, wherein said defined subset of CpG loci in said vimentin DNA comprises loci 37, 40, and
 45. 44. The method of claim 42, wherein said at least one marker DNA comprises BMP3 DNA, wherein said defined subset of CpG loci in said BMP3 DNA comprises loci 34, 53, and
 61. 45. The method of claim 42, wherein said at least one marker DNA comprises Septin9 DNA, wherein said defined subset of CpG loci in said Septin9 DNA comprises loci 59, 61, 68, and
 70. 46. The method of claim 42, wherein said at least one marker DNA comprises TFPI2 DNA, wherein said defined subset of CpG loci in said TFPI2 DNA comprises loci 55, 59, 63, and
 67. 47. The method of claim 42, wherein said at least one marker DNA comprises a EYA4 DNA, wherein said defined subset of CpG loci in said EYA4 DNA comprises loci 31, 34, 37, and
 44. 48. The method of claim 42, wherein at least one marker DNA comprises vimentin, BMP3, Septin9, and TFPI2 DNA. 